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I I We have obtained [Mg/Fe] measurements for 76.3% of the stars in the MILES 

P , spectral hbrary used for understanding stehar atmospheres and stehar populations in 

galaxies and star clusters. These abundance ratios were obtained through (1) a com- 
^ I pilation of values from the literature using abundances from high-resolution spectro- 

scopic studies and (2) a robust spectroscopic analysis using the MILES mid-resolution 
optical spectra. All the [Mg/Fe] values were carefully calibrated to a single uniform 
scale, by using an extensive control sample with results from high-resolution spectra. 
■ The small average uncertainties in the calibrated [Mg/Fe] values (respectively 0.09 and 

CO ' 0.12 dex with methods (1) and (2)) and the good coverage of the stars with [Mg/Fe] 

, over stellar atmospheric parameter space of the library will permit the building of 

new simple stellar populations (SSPs) with empirical a-enhancements. These will be 
available for a range of [Mg/Fe], including both sub-solar and super-solar values, and 
. for several metallicities and ages. These models will open up new prospects for testing 

' and applications of evolutionary stellar population synthesis. 
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1 INTRODUCTION history compared to the more quiet one of the solar neigh- 

^ , „ , , , „ bourhood (Tinsley 1980). 
Evolutionary stellar population synthesis, i.e. modelling 

spectral energy distributions emitted by evolving stellar An obvious solution, explored recently by several au- 

populations, is a natural approach to studying the stellar thors (e.g. Coelho et al. 2007; Walcher et al. 2009; Lee, 

content of different galaxies. One of the main ingredients of Worthey & Dotter 2009), would be to use theoretical li- 

these models are the stellar libraries, which can be empiri- braries with the desired coverage in chemical abundances, 

cal or theoretical. Stellar population models usually consider However, while theoretical libraries have improved dramat- 

only the total metal content of stars and, therefore, ignore ically in the last few years (Chavez, Malagnini & Morossi 

the different chemical abundance patterns that are present 1997; Murphy & Meiksin 2004; Rodriguez-Merino et al. 

in individual stars. However, different chemical abundance 2005; Munari et al. 2005; Martins et al. 2005; Coelho et al. 

patterns have a strong influence on the shape of the spectra. 2005; Fremaux et al. 2006; Bertone et al. 2008), they still do 

In particular, it is well known that stellar population models not reproduce real stars of all spectral types, with especial 

based on empirical libraries (which are mostly composed of problems depending on the wavelength range (i.e. Martins 

solar neighbourhood stars) cannot reproduce the high vab & Coelho 2007; Bertone et al. 2008). Some of the remaining 

ues of Mg abundances found in giant elliptical galaxies. This problems are the incompleteness of the atomic and molecu- 

is commonly interpreted as a consequence of high [Mg/Fe] lar hue opacity lists in the blue region of the spectrum and 

in these systems, most likely due to a rapid star formation for cool stars (T^ff < 4500 K) as well. 

Another approach is to compute, with the help of model 
atmospheres, response functions to characterise the varia- 

* E-mail: acmilone@das.inpe.br, andre.milone@inpe.br tion of specific spectral characteristic (usually Lick indices) 
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to variation of different elements (sec Trippico & Bell 1995; 
Korn, Maraston & Thomas 2005). Models using these re- 
sponse functions are those of Tantalo, Chiosi & Bressan 
(1998), Trager et al. (2000a), and Thomas, Maraston & Ben- 
der (2003). However, the accuracy of these theoretical pre- 
dictions have not been tested empirically yet. 

It is the main objective of this work to provide [Mg/Fe] 
abundance ratios for one of the most complete empirical 
stellar libraries currently available (MILES) (Mid-resolution 
Isaac Newton Telescope Library of Empirical Spectra, 
Sanchez-Blazquez et al. 2006). The MILES databeise, which 
was especially designed for stellar population modelling, 
contains flux calibrated optical spectra of high signal-to- 
noise ratio (S/N) for 985 stars covering AA3525-7500 A with 
a homogeneous resolution AA = FWHM = 2.3 A. The para- 
metric coverage of sample stars in the three-dimensional H- 
R diagram is quite wide: 2800 < T<,ff < 50400 K, 0.0 < log 
g < +5.0, and -2.7 < [Fe/H] < +1.0 dex, where [Fe/H] 
= log(Fe/H)* - log(Fe/H)0 such that formally log(Fe) = 
log(n(Fe)/n(H)) + 12, log(H) = 12 and n(Fe) and n(H) 
are the numerical densities (cm~^) of iron and hydrogen 
atoms respectively. The scales for these photospheric param- 
eters were carefully defined by Cenarro et al. (2007). Their 
precisions, respectively ±100 K, ±0.2 and ±0.1 dex, makes 
MILES good for SSP modelling. For the present work, we 
excluded those stars with uncertain or wrong atmospheric 
parameters (see Vazdekis et al. 2010 for details about how 
these stars were identified). 

The MILES [Mg/Fe] catalogue presented here consists 
of two measurement sets. The first one is composed of mear 
surements obtained from the literature from high spectral 
resolution analyses properly calibrated to a common system. 
The second set assembles abundances measured by us di- 
rectly from the MILES mid-resolution (hereafter MR) spec- 
tra and calibrated using the high-resolution (hereafter HR) 
sample. The catalogue is represented in two separate tables 
for field and cluster stars. Both tables are only available in 
electronic form. The paper layout is as follows: Section 2 
describes the compilation of HR abundance measurements 
from the literature and their calibration; Section 3 shows 
the Mg abundance measurements from the MILES spectra; 
Section 4 compiles the MILES [Mg/Fe] catalogue and anal- 
yses its coverage over the library parameter space; Section 
5 compares our stellar data with predictions of theoretical 
models focusing on the behaviour of some Lick indices with 
[Mg/Fe]; and, finally. Section 6 plans for applications to 
building new simple stellar population models with variable 
a-enhancement. Section 7 summarizes the whole paper and 
final conclusions. There are also three appendices: Appendix 
A confronts the compiled HR data with a well-known stellar 
spectrum library. Appendix B presents comparisons of the 
MILES photospheric parameter scales with those from the 
compiled HR studies, and Appendix C compares the results 
for cluster stars with HR studies. 



2 COMPILATION OF MAGNESIUM 

ABUNDANCES FROM HIGH-RESOLUTION 
STUDIES 

The first step of this work consisted of performing a biblio- 
graphic compilation of magnesium abundances from high- 



resolution spectroscopic analyses for the MILES library 
stars. To guarantee homogeneity between the measurements 
provided by several studies we performed a calibration and 
correction of systematic differences among sources and a 
chosen standard reference system, following a similar proce- 
dure as in Cenarro et al. (2001, 2007). For instance, Feltzing 
& Gustafsson (1998) give a detailed error analysis of elemen- 
tal abundances for G and K metal-rich dwarfs also including 
comparisons with other studies. In this section, we describe 
the chosen reference sample and the procedures we followed 
to homogenize the measurements to a single uniform scale 
of [Mg/Fe]. 

2.1 A reference scale for [Mg/Fe] 

Our reference sample to define a scale for the Mg/Fc abun- 
dance ratio is from Borkova & Marsakov (2005, hereafter 
BM05). The catalogue of Borkova & Marsakov is a robust 
compilation of the atmospheric parameters Tefi, log g and 
[Fe/H] plus [Mg/Fe] from high-S/N high-resolution analy- 
ses of field stars published between 1989 and 2003 (covering 
36 studies with Mg abundance determinations for around 
900 stars). BM05 computed weighted average values and 
their errors through an iterative procedure in order to cor- 
rect for the systematic deviations of each data set relative 
to reduced mean homogeneous scales. The uncertainty of 
[Mg/Fe] in BM05 is 0.05 and 0.07 dex, respectively, for 
metal-rich ([Fe/H] > —1.0 dex) and metal-poor stars ([Fe/H] 
^ —1.0 dex). The BM05 catalogue contains 218 stars in com- 
mon with the MILES library (all with log g ^ 3.0). 

BM05 was also the reference work adopted in Cenarro et 
al. (2009) , where Mg and Ca abundances were compiled for 
192 stars of their calcium triplet stellar library (hereafter 
CaT) of 706 objects. As there are many MILES stars in 
common with the CaT sample (132 stars), this work provides 
values that we can compare with (see Appendix A). 

2.2 Calibration of the high-resolution [Mg/Fe] 

We first checked the possible presence of systematic differ- 
ences in the scales of MILES and BM05 atmospheric param- 
eters ([Fe/H], log g and Tcff), however we did not find any 
(see Appendix B for details). Apart from the BM05 compi- 
lation, we obtained [Mg/Fe] for 97 more stars from 15 other 
HR studies, as listed in Table 1. Their abundance ratios 
were then carefully transformed onto the adopted scale as 
described next. 

Figure 1 shows the comparison of [Mg/Fe] for the stars 
in common between BM05 and other HR works. As can be 
seen, the relations are usually well described by an offset 
or a linear transformation. We derive these linear transfor- 
mations using a 3-cr clipping least-square {Isq) method min- 
imizing the distance in both axis (as the uncertainties in 
different studies are of comparable order). 

[Mg/Fc]work =A + B[Mg/Fe]BM05 (1) 

where [Mg/Fe]„ork represents the values computed in those 
HR works different from BM05. In BM05 and Cenarro et al. 
(2007), the comparison sample was gradually increased as 
each set of stellar parameters, like [Fe/H], was calibrated to 
a uniform scale. We, however, calibrate the [Mg/Fe] values 
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Table 1. Parameters of the linear calibrations of different sets of [Mg/Fe] to the BM05 scale and other information about the consulted 
high-resolution spectroscopic works. References are shown in the first column as: CGSOO (Caxretta, Gratton & Sneden 2000), FOO 
(Fulbright 2000), Ge03 (Gratton et al. 2003), BeOS (Besnby et al. 2005), T98 (Thcvcnin 1998), RLA06 (Reddy, Lambert & AUende 
Prieto 2006), LH05 (Luck & Heiter 2005), EN03 (Erspamer & North 2003), FG98 (Fcltziug & Gustafsson 1998), Ce02 (Caliskan et al. 
2002), FK99 (Fulbright & Kraft 1999), H02 (Heiter 2002), AcOl (Adclman et al. 2001), Ae06 (Adclman et al. 2006), and Cc09 (Ccnarro 
et al. 2009). The —A/B and l/B values (second and third columns) are, respectively, the additive and multiplicative coefficients of 
the calibration expressions (Eq. 2). When they are represented by integer numbers it means that no calibration of [Mg/Fe] was applied 
because there are too few stars or none in common between the work and BM05 samples to compute a linear fit. Nc (forth column) 
represents the number of stars in common between each work and BM05 after applying the excluding criterion to each Isq fit. N; (fifth 
column) is the number of stars from each reference that are eligible to be included into the MILES [Mg/Fe] catalogue. Nr (sixth column) 
gives the number of MILES stars that are repeated in other works (once in each case). The [Fe/H] range of the stars to be included 
into our catalogue is shown in the seventh column and their [Mg/Fe] ranges (from the original values) are shown in the ninth column. 
The [Fe/H] and [Mg/Fe] intervals of the work samples in common with, respectively, the MILES and BM05 catalogues can be read in 
the plots of Figs. 1 and Al. The typical uncertainties of [Fe/H] and [Mg/Fe] for each work are presented, respectively, in the eighth and 
tenth columns. In the eleventh column, the main physical constrains on measuring the Mg abundances are cited (LTE, non-LTE, and 
the ionization stage of the Mg lines used in the abundance determination). The propagated uncertainty of [Mg/Fe] over the calibration 
process for each data sample is written in the last column. The weighted averages of [Mg/Fe] errors are given in the eleventh row for 103 
stars whose data come from the first nine works. The weighted averages of [Mg/Fe] errors for 9 stars of the five last listed works (that 
do not have stars repeated in other works nor were their data calibrated to the [Mg/Fe] uniform scale in the current work) is shown in 
the penultimate row. 
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separately for each work, basically because the comparison 
sample adopted here is large enough (218 stars from the 
BM05 compilation) and because this avoids the error prop- 
agation through the transformations. The calibrated values 
are obtained, then, inverting the Eq. 1 as following: 

[Mg/Fe]HR = {-A/B) + (l/B)[Mg/Fe]work (2) 

The transformation was performed only when A and B were 
significantly difi'erent from and 1 respectively, based on the 
student t-test with a 95% confidence level. 

Table 1 presents the list of works and the number of 
stars in common between each sample and the BM05 cata- 
logue, as well as the number of MILES stars to be included 
into our catalogue from each HR study and the number of 



MILES stars of each work duplicated in other source(s), 
whose total is sixteen. The calibration constants and the 
typical work uncertainties of [Fe/H] and [Mg/Fo] are shown 
in this table. We also show the [Fe/H] and [Mg/Fe] ranges 
of the stars included from these works. Comments about 
the LTE assumptions and spectral lines employed in each 
work arc given in the table too. Statistically reliable linear 
calibrations of [Mg/Fe] were applied for the first nine works 
listed in Table 1. The data from Carretta, Gratton & Sneden 
(2000), Fulbright (2000), Bensby et al. (2005), Erspamer & 
North (2003), and Feltzing & Gustafsson (1998) wore not al- 
tered because there is no detectable difference between the 
[Mg/Fe] scales of these works and the BM05 scale. No cal- 
ibration was applied, either, to the data from the 6 last 
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Figure 1. [Mg/Fe]„ori; as a function of [Mg/Fe]BM05- 9 panels (a to i) showing comparisons between different samples with the reference 
set from Borkova & Marsakov (2005). The statistically representative, 3-cr-clipped (illustrated by parallel dashed blue lines), linear Isq 
fittings are presented by the thick blue lines (clipped data is represented by red symbols, occurring in the last panel only). The constants 
A and B and parameters rms and r (correlation coefficient) of the linear fits (Eq. 1) are listed at the top of each panel. The calibration 
expressions of [Mg/Fe],vork to [Mg/Fe]BM05 (Eq. 2) are also shown after applying the 95% t-tcst. The work designation is cited at the 
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entries of Table 1: Caliskan et al. (2002) and Fulbright & 
Kraft (1999), because there are very few stars in common 
with the BM05 sample; Heiter (2002) and Adelman et al. 
(2001, 2006) because there is no star in common; and Ce- 
narro et al. (2009) because their data are already on the 
same system we have adopted too (see Appendix A). 

Sixteen stars have [Mg/Fe] values from duplicated 
sources. When the difference between distinct sources is 
larger than 4a, then the values with smaller uncertainties 
were adopted (3 cases only). The final abundance ratios for 
the 13 remaining duplicated cases were computed as sim- 
ple averages after separately calibrating to the BM05 scale. 
These stars were also used to evaluate the calibration pro- 
cess as a whole and helped us to compare the uncertainties 
of [Mg/Fe] when there are duplicated or single data sources 
(see Sect. 2.3). 

Summarizing, the [Mg/Fe] catalogue of MILES contains 
315 stars with HR measurements (see Table 1): 

(i) 218 stars wfth [Mg/Fe] collected directly from BM05, 

(ii) 91 stars with [Mg/Fe] obtained from other published 
works and calibrated to the same single uniform scale, and 

(iii) 6 stars whose [Mg/Fe] ratios were collected from 
other works and inserted into it without any transforma- 
tion. 



2.3 Precision of calibrated [Mg/Fe] 

The uncertainties of calibrated [Mg/Fe] to the BM05 scale 
were estimated through the propagation of their original er- 
rors taking also into account the precision of calibration pa- 
rameters (see Table 1 and plots of Fig. 1). The error propa- 
gation through the linear expressions of calibration process 
was based on adding variances. The [Mg/Fe] uncertainties 
can be summarized as: 

(i) 0.07 and 0.05 dex for the 218 stars that define the base 
uniform scale, for respectively, [Fe/If] 5C —1.0 dex (47 stars) 
and [Fe/H] > -1.0 dex (171 stars) (as described in BM05); 

(ii) 0.17 dex for 75 stars (weighted average) whose abun- 
dance ratios were calibrated from single values; 

(iii) 0.10 dex for 13 stars based on average calibrated 
abundance ratios (data from duplicated sources); 

(iv) 0.14 dex for 3 stars from Ce09; and 

(v) 0.13 dex for 6 other stars whose [Mg/Fe] were not 
transformed onto the BM05 scale. 

The averaged [Mg/Fe] values compiled from duplicates 
are statistically more precise (la = 0.10 dex) than the cal- 
ibrated abundance ratios that have been computed from 
single sources {la — 0.17 dex). The final weighted aver- 
age uncertainty of [Mg/Fe] is around 0.09 dex over all com- 
piled HR data (315 stars). It is 0.10 dex for 88 metal-poor 
stars ([Fe/H] < -1.0 dex) and 0.08 dex for 227 metal-rich 
ones ([Fe/H] > — 1.0 dex). The mean uncertainty of [Mg/H] 
over whole range of metallicity is 0.13 dex, estimated by the 
quadratic sum of cr [Fe/H] and a [Mg/Fe]. 
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2.4 Extended control sample for calibrating 
[Mg/Fe] 

In total, we compiled [Mg/Fe] for 315 MILES stars covering 
about 1/3 of the library (263 dwarfs and 52 giants, respec- 
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Figure 2. [Mg/Fe] vs. [Fe/H] for the HR data of MILES [Mg/Fe] catalogue (summing 315 stars). [Fe/H] is on the MILES scale (Cenarro 
et al. 2007). The sources are represented by different symbols: BM05 catalogue as filled blue circles (218 stars), CGSOO as diagonal crosses 
(2 objects), FOO as four-pointed stars (13 objects), Ge03 as five-pointed stars (1 object), T98 as eight-pointed stars (35 objects), LH05 
as open squares (12 objects), EN03 as open pentagons (7 objects), FG98 as open hexagons (5 objects), FK99 as filled red triangles (1 
object), H02 as filled red squares (1 object), AeOl as filled red pentagons (3 objects), Ae06 as filled red hexagons (1 object), and Ce09 as 
filled black circles (3 objects). The data for 13 stars, whose sources are duplicated, are represented by open cyan circles. The red symbols 
represent stars whose [Mg/Fe] were not calibrated to the uniform scale (6 cases). The designations of works are cited in the caption of 
Table 1. Three chemically peculiar stars are identified by their names. 



lively 49% and 12% of them), which was used to define an 
extensive control sample for calibrating our own Mg abun- 
dance measurements at mid-resolution (Sect. 3). Dwarfs as 
designated when log g ^ 3.0 whilst giants when log g < 3.0, 
as in BM05 and the MILES database itself. 

Figure 2 shows [Mg/Fe] versus [Fe/H] for the HR part 
of the MILES [Mg/Fe] catalogue and its parametric cov- 
erage is shown in Table 6 (Sect. 4). Specifically, the con- 
trol sample owns 306 stars covering nearly 31% of MILES 
(255 dwarfs, and 51 giants) and practically presenting the 
same coverage. Besides six stars whose [Mg/Fe] has not been 
transformed onto the reference scale, three chemically pecu- 
liar metal-poor stars (BD-(-800245 confirmed by Ivans et al. 
2003, HD097855, and HD192640) have also been excluded 
from the control sample because their [Mg/Fe] lay far from 
the [Mg/Fe] vs. [Fe/H] trend described by the stars in the so- 
lar neighbourhood (see Fig. 2). Recently, Nissen & Schuster 
(2010) have classified a kind of galactic objects named low- 
a(-enhancement) stars distributed over two distinct nearby 
halo populations based on their kinematics, whose origins 
might due to the accretion from dwarf galaxies and some 



of them from the Omega Centauri globular cluster (denom- 
inated as a progenitor galaxy). HD097855 might belong to 
the proposed group. All of those nine stars have been incor- 
porated into our catalogue. 



3 MAGNESIUM ABUNDANCES MEASURED 
AT MID-RESOLUTION 

To extend the magnesium abundance characterization of 
the MILES stars, a spectroscopic analysis based on a LTE 
spectral synthesis of Mg features was carried out at mid- 
resolution. 

Recovering element abundances in stellar photospheres 
at medium spectral resolution (resolving power between 
1,000 and 10,000) has been a well-established and alterna- 
tive approach for many decades. For instance, Pagel (1970) 
reported that low-resolution spectroscopic analysis was one 
way to find the metallicity of nearby stars, by calibrating 
the results with a reference sample which, indeed, has been 
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highlighted as an important step involved (Friel & Janes 
1993; Kirby et al. 2009; Marsteller et al. 2009). 

Our present Mg abundances reported in the current sec- 
tion of the paper are based on the MILES spectra, which 
have a resolution comparable to many previous studies (e.g. 
Chavez, Malagnini & Morossi 1995; Terndrup, Sadler & Rich 
1995; Cook et al. 2007). 

3.1 Computation of the synthetic spectra 

Our analysis is based on a LTE spectral synthesis computed 
with the most recent stable version of the MOOG code at the 
time of developing this work (Sneden 2002). The synthesis 
code was fed by linearly interpolated model atmospheres 
over the MARCS 2008 grid (Gustafsson et al. 2008), up- 
to-date atomic line lists from the VALD database (Vienna 
Atomic Line Database, Kupka et al. 2000, Kupka et al. 1999, 
Ryabchikova et al. 1997, and Piskunov et al. 1995) and a set 
of important molecular lines of C2, CN and MgH (Kurucz 
1995) in order to compute a series of model spectra for each 
MILES star ranging over five values of [a/Fe]: —0.60, —0.30, 
0.00, -H0.30, +0.60 dex. 

The abundance of all Q-particle-capture elements - O, 
Ne, Mg, Si, S, Ar, Ca and Ti - were equally and simultane- 
ously modified to represent a global variation of their chem- 
istry in a stellar photosphere. Although these elements may 
have distinct individual abundance ratios relative to iron, 
assuming homogeneity for them should be a reliable approx- 
imation among many computations such as molecular disso- 
ciative equilibrium and partial pressure of different species, 
especially when the element abundances are not individually 
and previously known in each star. Moreover, changing glob- 
ally these ratios for all a-elements is coherent with the model 
atmospheres adopted in our spectral synthesis. We tested 
the alternative approach at medium spectral resolution (i.e. 
varying the Mg abundance only, keeping the abundances of 
other a?-elements fixed) for some MILES stars and the ef- 
fect was negligible. Even adopting the homogeneous global 
trend for the Q-elements, their individual abundances can be 
independently measured by separately analysing their own 
absorption features. For instance, the Mg abundance is only 
consistently quantified as one or more magnesium lines are 
analysed instead of computing and measuring features of 
other a-elements. 

We selected model atmospheres whose chemistry fol- 
lows the general pattern of solar neighbourhood for the a- 
elements (the standard composition group of the MARCS 
2008 models, Gustafsson et al. 2008), i.e. [a/Fe] = -f0.40 dex 
when [Fe/H] -1.00 dex and [a/Fe] =0.00 dex for [Fe/H] 
> 0.00 dex with intermediate variable values between these 
fiducial metallicities (-1-0.30, 4-0.20 and 4-0.10 dex for [Fe/H] 
= —0.75, —0.50 and —0.25 dex respectively). On this stan- 
dard range of compositions, every other element (X) follows 
the iron abundance, [X/H] = [Fe/H] meaning [X/Fe] = 
dex. A MARCS model presents a stratification over 56 plane- 
parallel layers for its photosphere. In the case of giants, the 
models (with -0.5 ^ log p ^ 3.5 in the MARCS grid) have 
been originally generated under a three-dimensional geom- 
etry for different masses and afterwards transformed onto 
one-dimensional representations by their own developers. 
We chose 3-d models with one solar mass for giants (log 
g < 3.0 in the current work). All chosen models for dwarfs 



(with 3.0 ^ log (7 ^ 5.5 in the grid) and giants have a micro- 
turbulence velocity of 2.0 km s"^ The MARCS 2008 grid 
ranges are: 

(i) 2500 < Tcft < 8000 K with steps of 100 K for TcH < 
4000 K and 250 K for Toff > 4000 K; 

(ii) —1.0 ^5 log g ^ 5.0 (or 5.5 in some cases) with 0.5 
constant steps; and 

(iii) [Fe/H] = -5.00, -4.00, -3.00, -2.00, -1.50, -1.00, 
-0.75, -0.50, -0.25, 0.00, +0.25, +0.50, +0.75 and +1.00 
dex. 

Linear interpolations of the model atmospheres were 
automatically done by using the user-friendly software 
of Masseron (2008), which is publicly available from the 
MARCS models' web-site ( |http://marcs. astro. uu.se/"] ). The 
layers of each input model to the MOOG code were repre- 
sented by the optical depth at 5000 A, the thermodynamic 
equilibrium temperature, the total gas pressure, and the 
electron numeric density (all quantities in cgs units). The 
micro-turbulence velocity was fixed at 2.0 km s~^ for all 
layers of each model. The effective broadening is dominated 
by the spectroscopic instrumentation and is suitably repre- 
sented by a Gaussian convolution. 

Our spectral synthesis computations adopt the same so- 
lar abundance pattern used to build up the MARCS 2008 
model atmospheres (Grevesse, Asplund & Sauval 2007, here- 
after GAS07), aiming at an important internal consistency 
regarding the reference solar chemistry. The solar abun- 
dances in GAS07 of those a-elements on a logarithm scale 
are: log(O) = 8.66, log(Ne) = 7.84, log(Mg) = 7.53, log(Si) 
= 7.51, log(S) = 7.14, log(Ar) = 6.18, log(Ca) = 6.31, and 
log(Ti) — 4.90. The abundance of iron is given by log(Fe) = 
7.45 on this scale. 

Many authors re-scale the oscillator strengths {gf) of 
atomic lines to reproduce the solar spectrum, i.e. cali- 
brate the line-strengths to the Sun's photospheric condi- 
tions. However, we preferred to maintain their laboratory gf 
values from VALD instead of normalizing them to the solar 
scale because we intend to calibrate our MR measurements 
using our HR control sample. The theoretical spectra were 
calculated at 0.02 A wavelength-steps by assuming an opac- 
ity contribution for the stellar continuum at 0.50 A bins. 
For the often strong Mg b lines (see Sect. 3.2 and Table 
2), we adopted the Unsold approximation for the interac- 
tion constant (Ce) of the van der Waals dump parameter 
76 multiplied by a 6.3 factor. This dumping represents col- 
lisions among neutral atomic species mainly H I and He I 
as in cold photospheric layers. On the other hand, the same 
constant multiplied by a MOOG internal factor was chosen 
for the other Mg feature. The Unsold approximation (de- 
fault in MOOG) means that Ce is basically due to the H I 
atoms as a function of the excitation potential of an elec- 
tronic transition, but 76 remains also dependent on the local 
gas pressure and temperature. 

All our model spectra were computed at the MILES 
resolution, which is slightly different in each Mg feature re- 
gion (FWHM = AA = 2.40 A for Mg5183 and 2.35 A for 
Mg5528), by applying Gaussian smoothing to represent the 
(instrumental dominated) broadening. At this resolution, no 
additional stellar rotational broadening needs to be consid- 
ered unless the line-of-sight rotational velocity Vi-otsin(i) is 
greater than ~130 km s~^. The wavelength scales of all 
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Table 2. The magnesium features: their central pass-bands 
Featurcband and two pseudo-continuum windows Cb and Cr, 
adopted for defining linear local continua for measuring their 
pseudo-equivalent widths. 



Feature 


Cb 


FeaturCband 


Cr 




(A) 


(A) 


(A) 


Mg5183 


5177.4-5178.3 


5179.2-5187.3 


5199.9-5200.8 


Mg5528 


5518.0-5521.0 


5524.0-5531.0 


5539.0-5542.0 



MILES observed spectra were carefully shifted to the rest 
wavelength to match the theoretical model scales, i.e. wave- 
lengths in air as available in the VALD database for A 
from 2000 A up to the infrared. Spectral cross-correlations 
were applied for this purpose through the cross-correlation 
fxcor task of the Radial Velocity Analysis Package of the 
NOAO Optical Astronomy Packages of IRAF (Image Re- 
duction and Analysis Facility |3 by adopting a given corre- 
spondent model spectrum for each star as a template. The 
model spectra computed for either [Mg/Fe] — -1-0.60 dex or 
-1-0.30 dex were chosen in order to guarantee reliable spec- 
troscopic cross-correlations for all MILES spectra, from the 
hotest metal-poor to the coldest metal-rich ones, i.e. with 
strong absorption lines in the template spectra. Then the 
wavelength scales of all computed spectra were re-binned 
to exactly agree with the sampling of the MILES spectra at 
each observed wavelength bin (0.90 A). Finally, flux normal- 
izations of the empirical spectra based on the local pseudo- 
continuum were carefully applied at the regions around each 
Mg feature in order to adequately match to the flux normal- 
ized scale of the theoretical spectra. Fiorentin et al. (2007) 
emphasized the importance of performing reliable compar- 
isons on compatible flux and wavelength scales between ob- 
served MR spectra and synthetic ones, paying special atten- 
tion to choosing useful absorption features and accurately 
calibrating the abundance results. 



3.2 Mg features analysed: Mg5183 and Mg5528 

Two strong Mg features were chosen and carefully tested to 
be measurable and useful for recovering magnesium abun- 
dances at the MILES spectral resolution with acceptable 
precision for generating new stellar population models with 
[Mg/Fe] constraints, i.e. with comparable precision to that 
of [Fe/H] in MILES (0.10 dex). These features are: (i) the 
reddest line of the Mg b triplet (A5183.604 A) named here 
Mg5183, which is usually the strongest and the most Mg- 
sensitive of the three lines, and (ii) MgIA5528.405 A, here- 
after Mg5528. 



^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation, USA. 



3.3 Methods applied: pseudo-equivalent width 
and line profile fit 

Two methods were chosen to measure the magnesium abun- 
dances: 

(i) based on pseudo-equivalent widths, hereafter EW 
(pseudo-ones in fact due to the extensive line blanketing 
at the MILES spectrum resolution), and 

(ii) applying line profile fittings, hereafter LPF. 

These methods are usually adopted on HR and MR 
analyses and both require a previous knowledge of the pho- 
tospheric parameters. When there are more than 2 lines of an 
element, the modelling of their equivalent widths are prefer- 
able instead of fitting their profiles. The fit of line profiles 
is commonly adopted to extract element abundances from 
molecular absorptions where a myriad of lines from a sin- 
gle substance are very close to each other. When analysing 
atomic features, both methods generally need isolated lines, 
however, they can work on composite lines that can be de- 
blended in single profiles if the abundances of the others 
absorbers are already known. Traditionally the EW method 
is applied in an automatic process to extensive sets of lines. 
On the other hand, the second method is employed focusing 
on careful visual inspection in a feature-by-feature and star- 
by-star base. The abundance precision derived from each 
method depends on how many features are adopted and 
how sensitive each absorption line is to a given elemental 
abundance variation. 

At each chosen Mg feature region, we automatically ap- 
plied these methods as explained below. 

(i) The equivalent width of a weak line is propor- 
tional to the number of absorbers of element X (EW cx 
n(X)), whilst the equivalent widths of strong lines are de- 
pendent on n(X)^^'^. Therefore direct comparisons were 
made on planes [Mg/Fe] vs. log(EW) by adopting sim- 
ple linear Isq fittings for each set of theoretical equiv- 
alent widths. Note that the iron abundance [Fe/H] is 
fixed in all model computations for each star (like the 
other photospheric parameters assuming those compiled by 
Cenarro et al. 2007) and [Mg/Fe] = [Mg/H] - [Fe/H]. 
Therefore we adopt that relationship instead of the direct 
log(EW) vs. log(n(X)). Consequently, [Mg/Fe] can be di- 
rectly measured instead of [Mg/H] and they were obtained 
through interpolation of the [Mg/Fe] vs. log(EW) relation- 
ship for each star and feature combination. The equiva- 
lent width measurements of both Mg features were per- 
formed with the LECTOR code (A. Vazdekis' webpage, 
www.iac.es/galeria/vazdekis/SOFTWARE/) and INDEXF 
software (www.ucm.es/info/Astrof/software/indexf), how- 
ever, the uncertainties have been estimated through IN- 
DEXF only. Both software provided the same results within 
an accuracy of milli- Angstroms. The measurements were 
carried out within the central passband of each Mg feature 
by adopting a linear local flux continuum that is deflned by 
the average fluxes and wavelengths of two pseudo-continuum 
windows held very near to the passband (at each side of 
the feature). The central passband and pseudo-continuum 
windows of both Mg features were carefully chosen in or- 
der to provide representative measurements for their equiv- 
alent widths and are presented in Table 2. We computed 
the equivalent width uncertainties as being dominated by 
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Figure 3. Examples of spectral synthesis of the Mg5183 and Mg5528 features for a MILES' dwarf (HD006834: T^ff = 6473 K, log 
g = 4.22, [Fe/H] = -0.89, at the two top sub-figures) and giant (BD+182890: Toff = 4957 K, log g = 2.20, [Fe/H] = -1.61, at the 
two bottom sub-figures). Their central bandpass and two pseudo-continuum windows, Cj, and c^, are shown in the main panel of each 
sub-figure. The observed spectrum is represented by the solid black line and the synthetic ones by the colour dashed lines. Each set 
of theoretical spectra was computed for five a-enhancements ([a/Fe] = —0.60, —0.30, 0.00, -1-0.30, +0.60 dex, respectively, on colours 
blue, green, yellow, magenta and red) assuming the star's photospheric conditions fixed. Two small graphs illustrate the two abundance 
determination methods in each sub-figure: EW (bottom left corner) and LPF (bottom right corner) presenting their [Mg/Fe] 
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photon statistics based on the formahsm from Cardiel et al. 
(1998). The relative uncertainties are distributed between 
few percents and 40%, whilst the equivalent widths range 
from 0.15 up to 2.60 A and from 0.05 up to 0.95 A for the 
Mg5183 and Mg5528 features respectively. Wc noted that 
the EW uncertainties are dependent on the spectral S/N, 
as expected, and it is more evident for the Mg5183 feature. 
The higher S/N, smaller the EW relative error is. The S/N 
has been computed as an average between the blue and red 
regions of MILES stellar spectra, and it ranges from 10 up 
to 550 per A with typical values around 235 per A. It IS 
also noted there is a strong inverse correlation between the 
relative errors and the equivalent width values. 

(ii) The line profile comparison between the observed 
Mg feature and each corresponding set of synthesized fear 
tures was made within the central passband through rms 
statistics. Optimal [Mg/Fe] from each feature were derived 
for each MILES' star through rms minimization, as is 
shown in Fig. 3 for a dwarf and giant spectra. The pseudo- 
continuum windows were adopted to accurately match the 
continuum fluxes of the observed and model spectra. Careful 
attention was paid to specific cases for which the theoretical 
and observed spectra do not match each other simultane- 
ously at both continuum windows. Minor corrections were 
applied to the observed spectrum continuum flux for this 
purpose based on eye-trained inspections (additive or mul- 
tiplicative corrections). Each curve of rms as a function of 
[Mg/Fc] was fitted by a spline. 

Figure 3 shows examples of spectral synthesis of the 
Mg5183 and Mg5528 regions for a dwarf and giant. The 
main panels of each sub-figure (for each combination star- 
feature) present the MILES spectrum compared with model 
spectra. The sub-figures also include graphs illustrating the 
abundance measurement methods. 

The global procedure for obtaining and calibrating 
[Mg/Fe] ratios from our MR spectra is as follows: 

(i) automatic measurements at each Mg feature by apply- 
ing both methods to obtain [Mg/Fe] 

(ii) for each feature, computation of simple averages from 
the two methods when both show reliable results checked 
through visual inspections to get [Mg/Fe] feature (in some 
cases, only one method results in a reliable value and in 
other cases, both methods fail); 

(iii) linear calibration of [Mg/Fe] feature to a uniform scale, 
separately for each feature, via comparison with the HR 
star control sample (described at Sect. 2.4), to compute 
[Mg/Fc]F,1i'i,,; 

(iv) and simple averaging of the calibrated ratios obtained 
from both features combined when possible, in order to ob- 
tain the final calibrated abundance ratios [Mg/Fe]'^^''''. 

When the spectral synthesis of a given Mg feature does 
not satisfactorily reproduce the observed stellar spectrum 
(equivalent width and/or absorption line profile), wo call 
this a non-reproduced spectral synthesis and classify it such 
as (a) inadequate reproduction of the observed spectral con- 
tinuum, (b) saturation efi'ect on the [Mg/Fe] vs. log(EW) re- 
lationship (non-linearity), (c) extrapolation on the [Mg/Fe] 
vs. log(EW) relationship and/or on the rmsLPF vs. [Mg/Fe] 
one, (d) possible inaccurate model atmosphere interpolation 
(i.e. done around the borders of the MARCS 2008 grid), 
(e) possible higher rotation velocity than 130 km s~^, (f) 



absence of molecular lines to compute the model spectrum 
(e.g. TiO bands), (g) low quality observed spectrum (S/N 
below 50 per A), (h) possible wrong photospheric parame- 
ters, (i) suspect chemically peculiar star, and (j) other un- 
known causes and effects. Consequently, the abundance ratio 
provided by the correspondent method is not reliable. 

Calibration of each [Mg/Fe] feature WaS applied after av- 
eraging the measurements obtained from both methods. We 
decided to compute simple averages because each method 
explores one aspect of a reliable spectral synthesis. Whilst 
the EW comparison focuses on the reproduction of the total 
energy absorbed by the feature relative to the local contin- 
uum, the line profile fit takes into account the line shape 
including the core and wings. 

Moreover the absolute differences between EW and LPF 
methods are always smaller than 3cr [Mg/Fe] feature (around 
0.03 dex for Mg5183 and 0.07 dex for Mg5528). We checked 
for dependences of the differences between methods on the 
atmospheric parameters Tefi, log g and [Fe/H] and the re- 
sults do not show any dependence, for both features. We 
also noted there is no parametric dependence of the differ- 
ences ( [Mg/Fe] fellure - [Mg/Fe]HR) and ( [Mg/Fe] J^^f^.e - 
[Mg/Fe]HR) for both Mg features and methods. 



3.4 Calibration of the mid-resolution [Mg/Fe] 

The calibration to the previously adopted uniform scale re- 
lies on an inverse linear transformation, 

[Mg/Fe]?e1lure = {-a/h) + (1/6) [Mg/Fe] feature (3) 

that is obtained through a 3-(j clipping simple linear Isq fit 
of our MR measurements as a function of the HR calibrated 
values from the [Mg/Fe]HR star control sample (Sect. 2.4). 
We also applied a 95% t-test to verify if the fit parameters a 
and b are distinguishable from zero and unity respectively. 
In case they were not, wc would adopt a = and 6=1. 

Figure 4 shows the comparisons of [Mg/Fe] Mgsiss and 
[Mg/Fe]MgBB28 against [Mg/Fe]HR, in which the linear Isq 
fittings are presented simultaneously for dwarfs and giants. 
The computed calibration expressions are: 

[Mg/Fe]Sg5i83 = 0.058 dex + 1.546[Mg/Fe]Mg5i83 (4) 

[Mg/Fe] Sg^;.28 = 0.227 dex + 1.011[Mg/Fe]Mg5528 (5) 

The cahbrated abundance ratios [Mg/Fe]'''''''' from the 
Mg5183 and Mg5528 features that respectively lay out- 
side the intervals [-0.42 dex, +0.92 dex] and [-0.41 dex, 
-1-0.88 dex] are based on extrapolations over the calibra- 
tions. The intervals were estimated taking into account the 
[Mg/Fe]feature scale coverages of the stellar common sample 
between the HR compilation data and our MR measure- 
ments as well as the internal uncertainties of [Mg/Fo]Mg5i83 
and [Mg/Fe]Mg5528 (described in Sect. 3.5). Therefore the 
abundance ratios become more uncertain outside those in- 
tervals (see Fig. 4). 

The stars' samples adopted in the calibrations are ex- 
tensive enough and exhibit wide ranges of photospheric pa- 
rameters, which are as extensive as those of the MILES 
control sample, including uncertainties. Parameter intervals 
covered can be seen in plots of Fig. 5 (described ahead) . 
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Figure 4. [Mg/Fe]Mg5i83 vs. [Mg/Fe]HR, top panel (a), and 
[Mg/Fe] Mg5528 vs. [Mg/FeJiiR, bottom panel (b), comparisons 
for dwarfs and giants together in each, respectively represented 
by black filled circles and red filled triangles. The straight line 
inversely derived from the simple linear Isq fit [Mg/Fe] feature = o. 
+ 6[Mg/Fe]£fii is also shown for each comparison (blue solid line) 
with two parallel blue dashed lines illustrating the 3-cr clipping 
procedure. The parameters a and b are displayed on the top of 
graphs as well as the calibration expressions themselves. 



We verified the atmospheric parameter sensitivity of the 
[Mg/Fe] differences between MR and HR results. Figure 5 
presents them as a function of TcS, log g and [Fe/H] for 
both Mg features. There is no significant dependence on the 
photospheric parameters, as indicated by the hnear fits (tiny 
linear correlation coefficients). Furthermore the differences 
are always smaller or comparable to rms scatter along the 
whole parameter scales. 



Therefore, where possible, we averaged the calibrated 
ratios obtained from both features, in order to obtain cali- 
brated abundance ratios as follows: 

[Mg/Fe]-"b = ([Mg/Fe]S^5^83 + [Mg/Fe]^Mt"528)/2 (6) 

In Fig. 6, the residuals between calibrated results from 
the two features are plotted as a function of T^s, log g and 
[Fe/H]. There is no systematic dependency on the stellar 
parameters, as it is noticed from the small linear correlation 
coefficients r. This supports our procedure of computing av- 
erages of [Mg/Fe] jjj^turc- The medium data dispersion (0.15 
dex rms) is comparable with the typical internal error and 
systematic error for the calibrated ratio from Mg5528, but 
is slightly greater for the other feature. It is equal to the sys- 
tematic error of [Mg/Fe]Mg5528 (±0.15 dex), but it is slightly 
greater than cr[Mg/Fe]Mg5''i83 (0-13 dex) and cr [Mg/Fe] ^^'"b 
(0.10 dex). 

The criterion for deciding if the calibrated [Mg/Fe] 
measurements of both Mg features can be averaged to ob- 
tain representative determinations is fixed by the maxi- 
mum deviation acceptable between them, i.e. the condi- 
ti on is express ed by [[Mg/Fej^^^g^sag - [Mg/Fej^^s^ga | < 

4cr [Mg/Fe] '^''"b or A[Mg/Fe]A/95528-M<,5i83 ^ 0.40 dex, as 
illustrated in Fig. 6 and adopted for calibrated HR mea- 
surements collected from duplicated sources (Sect. 2.2). See 
Sect. 3.5 to get details about how (7[Mg/Fe]'==''''^ was esti- 
mated. The Mg5183 determination was adopted for the stars 
that do not follow this condition (5 dwarfs and 10 giants), 
because this feature provides a better precision. 

After separately calibrating [Mg/Fe] Mg5i83 and 
[Mg/Fe] Mg5528, using the extensive control sample, we were 
able to compute average or individual feature values for a 
great number of dwarfs and giants. In total we measured 
[Mg/Fe] at MR for 437 extra MILES stars (150 dwarfs 
and 287 giants). This represents around 44% of the whole 
spectral library. 



3.5 Uncertainty of the mid-resolution 
measurements 

The internal uncertainties of [Mg/Fe] are due to the MILES 
photospheric parameter errors and the abundance determi- 
nation methods applied (EW and LPF). 

We estimated the error propagation from the impreci- 
sion of Tcft , log g and [Fe/H] by adopting model atmospheres 
directly collected from the MARCS 2008 grid, whose pa- 
rameters coincide with the photospheric parameters of all 
MILES stars possible considering the typical errors (Icr). 
This sample of MILES stars is composed of 135 dwarfs and 
31 giants, and it has an excellent coverage in the MARCS 
model parameter space (distributing between 4200 and 6200 
K in Toff, 1.0 and 4.7 in log g and -2.10 and -1-0.40 dex in 
[Fe/H]). Then, we computed their abundance ratios based 
only on the EW method by using each Mg feature, here- 
after simply [Mg/Fe] atm- Afterwards, we compared them 
with those measurements obtained by adopting interpolated 
model atmospheres, whose parameters actually match Toff^, 
log g and [Fe/H] of the stars, naming them [Mg/Fe]interpoi. 
In this procedure, all [Mg/Fe] are not calibrated to the HR 
uniform scale, since only differences are needed for estimat- 
ing errors. The deviation ([Mg/Fe]atm — [Mg/Fe] Interpol) is 



12 A. de C. Milone, A. E. Sansom and P. Sdnchez-Bldzquez 



1.0 
0.5 



fc^ 0.0 

DO 

^-0.5 



A[Mg/Fe]„g.„„=(3.1Ze-05±1.16e-05)T,„-0.2B2±0.086 (r = 0.1710, rnis = 0.105) 



-1.0' 

4000 



l.Or 



6000 



- A[Mg/Fe]„„_„„-(0.031±0.008)log g-0.a25±0.030 (r=0.2547, rms=0.103) 



K 0.5 - 

3 

<u 

^ 0.0 - 
QO 

< 

-1.0 - 


1.0 

^ 0.5 

0.0 




2.0 3.0 
log g 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 

a[Mg/Fe]^R_„„-(0.053±0.009)[Fe/H]-0.077±0.008 (r-0.3484 


, rms-O.lOO) 




Mg5183 



-1.0 

[Fe/H] 



A[Mg/Fe]„g.„„=(-1.17e-06±1.59e-D5)T^„-0.223±0,069 (r= -0.0045. rms=0.151 



8000 4000 



l.Or 



6000 
Te„ (°K) 



- A[Mg/Fe]„„_„„-(-0.016±0.010)log g-0. 1 69±0.040 (r=-0.093a, rms=0.150) 



K 0.5 - 

3 

0) 

^ 0.0 -- 

< : I 

-i.ot — ^ 

0.0 



Mg552B 



1.0 : 

i 0.5 : 

^ 0.0 - 



A[Mg/Fe]„g_„„-(-0.052±0.012)[Fe/H]-0.360±0.01 1 (r--0.a579, rms-0.146) 



Mg5528 



-1.0 

[Fe/H] 



Figure 5. [Mg/Fe] differences between non-calibrated mid-resolution (MR) and calibrated high-resolution (HR) measurements 
(A[Mg/Fe]MR-HR = [Mg/Fe] feature - [Mg/Fe]HH,) computed for the features Mg5183 (left panel) and Mg5528 (right panel) as 
a function of Tefi, log g and [Fe/H] distinguishing dwarfs (black filled circles) and giants (red filled triangles). Simple linear Isq fittings 
between the differences and each photospheric parameter are also shown (blue solid lines). The resultant expressions, rms and correlation 
coefficient of these linear fits are displayed on the top of each graph. 



well distributed around zero over all parameter scales and 
it does not show any stellar parameter dependence. These 
differences at la level gave us a good estimation for the 
global internal uncertainty due to the star's parameter er- 
rors, (T[Mg/Fe]atm, which is around 0.09 dex for the Mg5183 
feature and 0.10 dex for Mg5528. 

The internal uncertainties in [Mg/Fe] due to the EW 
method itself, cr [Mg/Fe] ew, for both Mg features were esti- 
mated by computing the propagation of the observed equiv- 
alent width errors through the relationship [Mg/Fe] vs. 
log(EW). We noted that the EW uncertainties dominate 
over the propagation error from the [Mg/Fe] vs. log(EW) 
theoretical fits. The median of the ct [Mg/Fe] ew asymmet- 
ric distribution for each Mg feature characterises the typical 
uncertainty that it is around 0.09 dex for Mg5183 and 0.20 
dex for Mg5528. 

We also investigated the dependence of cr [Mg/Fe] ew on 
the atmospheric parameters. There is some correlation with 
[Fe/H] for Mg5183 only, where cr[Mg/Fe]EW is higher on 
average for metal-poor stars (~0.2 dex) in comparison with 
metal-rich stars (~0.1 dex) with a transition limit around 
[Fe/H] = —1.0 dex. Although a metallicity dependence is 
also noticed in HR measurements, as in BM05, we decided 
not to estimate a [Mg/Fe] ew as a function of [Fe/H]. 



The uncertainty from the LPF method depends on how 
accurately the minimum of the curve is determined (see Fig. 
3). Typically this is estimated to be between 0.05 and 0.10 
dex, giving a mean of cr[Mg/Fe]LPF around 0.075 dex. The 
final uncertainty due to the measurement process was com- 
puted as mean value of the abundance uncertainties ob- 
tained with the two methods, because, although they are 
not completely independent, the data is treated in differ- 
ent ways through them. Typical value of (T[Mg/Fe]j^^^j^^^ is 
0.06 dex for Mg5183 and 0.11 dex for Mg5528. The typical 
final internal error for each Mg feature, 0.11 and 0.15 dex 
for Mg5183 and Mg5528 respectively, was obtained as the 
quadratic sum of cr [Mg/Fe] atm and o" [Mg/Fe] ^^^^^j^^^. 

The systematic uncertainties of calibrated MR [Mg/Fe] 
were estimated by comparing these values directly with HR 
values, such that the rms of deviation ([Mg/Fe] feature ~ 
[Mg/Fe]HR) represent a good estimation, i.e. for each Mg 
feature we have 

iV 

a[Mg/Fe]^^"'' = (^ ^([Mg/Fe]^^'"' - [Mg/Fe]HR)')^''' (7) 
1=1 

where means the number of stars in the MR vs. HR 
comparisons for each Mg feature determination. Therefore 
[Mg/Fe] recovered by Mg5183 holds an uncertainty of 0.13 
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Figure 6. [Mg/Fe] differences between tiie cali- 

brated measurements obtained witii tiie two Mg features 
(A[Mg/Fe]Mg5528-Mg5i83 = [Mg/FeJ-^b^g - [Mg/Fe]-«\33) 
as a function of Toff, log g and [Fe/H] for dwarfs (black filled 
circles) and giants (red filled triangles). Least-square linear 
fittings between the differences and each photospheric parameter 
are also shown (solid blue lines) considering the data shown 
as filled symbols. Tw o parallel do tted black lines are drawn 
to representing the Acr[Mg/Fe]'^'^^^^ as a maximum acceptable 
deviation for A[Mg/Fe]Mg5528-Mg5i83 (details in the end of 
Sect. 3.5). The black open circles represent the dwarfs and red 
open triangles the giants for which average abundance ratios 
were not computed based on the 4(t criterion. These stars are 
not considered in the Isq fittings. 



4-0.40) with the alpha-poor ([a/Fe] = 0.00) and standard 
{[a/Fe] = -1-0.40) models for this metallicity, and (iii) (0.00, 
0.00) and (iv) (0.00, -1-0.40) with the standard {[a/Fe] = 
0.00) and alpha-enhanced ([a/Fe] — 4-0.40) models for this 
metallicity. The syntheses were applied to four stellar evolu- 
tion stages: main sequence at Tcfi = 5000 K and log g — 4.5, 
turn-off main sequence at TcH — 6000 K and log g — 4.0, 
sub-giant at Teg = 5000 K and log g — 3.5, red giant branch 
at Toff — 4500 K and log g = 1.0. However, just the red gi- 
ant stage could be adopted for the solar metallicity due to 
a limitation of the MARCS alpha-enhanced class coverage. 
Afterwards, we directly compared the equivalent widths of 
each Mg feature measured on theoretical spectra that have 
been computed for each metallicity-a-enhancement combi- 
nation but assuming two different model classes. 

At [Fe/H] = —1.50 dex, it is noticed that the EW varia- 
tions of both Mg features are very acceptable within typical 
uncertainties, i.e. < 8% for the main sequence and sub-giant 
stages and < 3% for the turn-off and red giant stages that 
are translated into [Mg/Fe] abundance ratio changes smaller 
than around 0.10 and 0.05 dex respectively. At solar metal- 
licity, the EW variations of red giant stage (the only one 
analysed) are smaller than 7% for Mg5183 (or about 0.09 
dex in [Mg/Fe]) and 20% for the Mg5528 feature (or ~0.2 
dex in [Mg/Fe]). Therefore it is viable to perform spectral 
synthesis at MR by fixing the model atmosphere chemistry 
and changing the Q-element abundances to cover a large 
range of [a/Fe] values as it was done in the current work. 

However, a conservative variation of ±0.2 dex for [a/Fe] 
in the spectral syntheses done under fixed model atmo- 
sphere chemistries could be assumed in order to figure out 
which abundance ratio determinations would be considered 
as extrapolations based on the a-enhancement compatibil- 
ity. Concerning the MARCS standard models adopted and 
taking into account the uncertainties of [Mg/Fe] , we have 
found 33 stars, which represent just 7.5% of all MR deter- 
minations and are identified in the catalogue (Sect. 4). Con- 
sequently, the [Mg/Fe] of each case might have a less accu- 
racy in the sense of it was based on a model atmosphere 
whose chemistry does not exactly follow the abundances of 
Q-elements adopted in the spectral synthesis. These cases in- 
clude those extrapolations over the Mg feature calibrations 
themselves (Sect. 3.4). 



dex or 0.15 dex when is uniquely measured by Mg5528. 
When it was possible to compute [Mg/Fe] as an aver- 
age from both feature determinations, its systematic error 
a[Mg/Fe]'^^'''' was estimated by quadratic mean reaching 
0.10 dex. Table 3 summarizes the [Mg/Fe] uncertainties. 

We also investigated the influence of the models atmo- 
sphere alpha-enhancement chemistry on the spectral synthe- 
sis carried out at MR by always keeping the model chem- 
istry unchanged at each metallicity (see Sect. 3.1). Besides 
the MARCS standard models adopted in our work, there 
are other model classes: one named alpha-poor assuming 
a/Fe solar for -2.00 ^ [Fe/H] -0.25 dex, and another 
denominated alpha-enhanced with [a/Fe] = -1-0.40 dex cov- 
ering [Fe/H] from —0.75 to -1-0.50 dex. Thus we performed 
spectral syntheses to cover four ([Fe/H], [a/Fe]) (dex, dex) 
combinations as following: (i) ( — 1.50, 0.00) and (ii) ( — 1.50, 



3.6 Coverage of the MR measurements 

There are 843 MILES stars with catalogued atmospheric 
parameters within the MARCS 2008 grid. From those, 308 
already had [Mg/Fe] measurements from HR studies (Sect. 
2). The 535 residual stars were spectroscopically analysed 
by us at MR. The MR determinations cover wide ranges in 
atmospheric parameters (see Figs. 8 to 11 described in Sect. 
4), reaching 81.7% efficiency or completeness level inside the 
MARCS parameter space when the stars with HR data are 
not considered. Rather than be uniform over all scales, our 
measurements actually complement the HR data. However, 
depending on the region, one Mg feature works better than 
the other. h Section 4 presents a more extensive discussion 
about the parametric coverage of all MR and HR determi- 
nations. 
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Table 3. Typical uncertainties of [Mg/Fe] (dex unity) from our mid-resolution measurements, showing the internal errors for the Mg5183 
and Mg5528 features due to the photospheric methods adopted (pseudo-equivalent widths, EW, and line profile fitting, LPF) and stellar 
atmospheric parameter errors. The total uncertainty of non-calibrated [Mg/Fe] feature ^■iid the systematic error of calibrated [Mg/Fe] 
for each Mg feature are presented respectively in the last and penultimate rows. 



# 




Mg5183 


Mg5528 


both 


Notes 






(dex) 


(dex) 


(dex) 




1 

2 


(T[Mg/Fe]EW 
o-[Mg/Fe]LPF 


0.09 
0.075 


0.20 
0.075 




from the EW method 
from the LPF method 


3 


(T[Mg/Fe]inethod 


0.06 


0.11 




from rows 1 and 2 variance averaged 


4 


fT[Mg/Fe]atm 


0.09 


0.10 




due to the photospheric parameter errors 


5 


0-[Mg/Fe] feature 


0.11 


0.15 




internal errors from rows 3 and 4 combined in quadrature 


6 


a[Mg/Fe]-iib^^ 


0.13 


0.15 


0.10 


systematic errors (from HR comparisons), variance averaged in the last column 



In general, the Mg5183 feature does not work well on the 
coldest stars (Tcff < 4000 K) due to the presence of strong 
molecular absorptions of TiO and MgH as well. Mg5528 is 
not satisfactorily applied for the hottest giants (Teff > 5500 
K) , since this feature becomes too weak and it is practically 
insensitive to abundance variation in these cases. In par- 
ticular, Mg5183 gives reliable abundance measurements for 
dwarfs and giants with temperatures between 4000 and 8000 
K, while Mg5528 basically works on dwarfs with 3500 up to 
8000 K. However, for giants, Mg5528 can only be applied 
with great confidence for 3600 < Teff ^ 5500 K. Both Mg 
features can be used in the whole metallicity range. 

In specified cases, the spectral synthesis does not work 
on a single or both features, due to line saturation, non- 
reproduction of spectrum continuum, extrapolation on the 
[Mg/Fe] vs. log (EW) and/or rmsLPP vs. [Mg/Fe] rela- 
tionship, incompleteness of line lists (mainly TiO bands), 
and low quality spectra in some cases. The continuum and 
EW/LPF extrapolation cases do not exhibit any stellar 
parametric dependence for Mg5183. On the other hand, we 
noted the cases of line saturation and incomplete line list 
occur in metal-rich cold stars ([Fe/H] > —1.0 dex with Tefi 
< 4000 K) for both features. A few of those non-reproduced 
spectral cases (that sum in total 79 dwarfs and 115 giants for 
Mg5183 plus 60 dwarfs and 97 giants for Mg5528) could be 
fixed, excepting those due to the incompleteness of line lists: 
8 cases for Mg5183 (5 dwarfs and 3 giants), and 15 cases for 
Mg5528 (5 dwarfs and 10 giants). The line saturation effect 
was solved by adopting a spline fit on the EW method and 
sometimes a smaller number of models (9 cases in total) . 



4 THE MILES MAGNESIUM ABUNDANCE 
CATALOGUE 

We have obtained [Mg/Fe] covering a bit more than 3/4 
(more precisely 752 stars or 76.3%) of the MILES stellar 
spectrum library (411 dwarfs and 341 giants, respectively 
around 76% and 77% of their totals) that arc suitable for 
SSP modelling, i.e. the typical systematic uncertainty of 
[Mg/Fe] is 0.105 dex on average over our whole catalogue. 
The stars' coverage in the four-dimensional parameter space 



of MILES Tcfi , log g, [Fe/H] and [Mg/Fe] is extensive, as dis- 
cussed in this section. If we only consider the MILES stars 
with complete sets of photospheric parameters that sum 946 
objects (Ccnarro et al. 2007), the coverage reaches 79.5%. 

The compiled [Mg/Fc] catalogue of MILES is presented 
in Tables 4 and 5 for field and cluster stars respectively. 
The whole catalogue is only published in electronic form. 
Basically, the catalogue tables provide the [Mg/Fe] values 
with their individual errors together with the sources from 
where they have been obtained (i.e. the reference in case 
it has been compiled from HR works or the Mg feature(s) 
when it is a MR measurement). 

Figure 7- (a) shows the stars' coverage over the modi- 
fied H-R diagram log g vs. Teff, Figure 7-(b) explores the 
coverage around the projection [Fe/H] vs. Tcef, and Figure 
7-(c) presents it on the plane log g vs. [Fe/H]. The MILES 
stars with and without [Mg/Fc] arc represented in all these 
plots, in which the parametric extension of MARCS 2008 
grid is also drawn. The stars with [Mg/Fe] are distinguished 
according to the origin of their measurements. Figure 8 ex- 
hibits the distribution of MILES [Mg/Fe] catalogue over the 
Teff, log g and [Fe/H] scales. 

The sum of all measurements (HR plus MR) show fairly 
flat distributions across the photospheric parameter space 
covered by the MARCS models (histograms of Fig. 8). The 
[Mg/Fe] HR data are distributed along the main sequence 
(MS) basically from Tea about 4500 up to around 10000 K 
and on the giant branch mainly from 4000 up to 5500 K, 
as seen in Fig. 7- (a). Our MR measurements have a wide 
distribution over the plane log g vs. Teft, however there are 
some deficiencies such as in the low-MS, the red giant branch 
tip and the hottest giants. The histogram of Tes , Fig. 8-(a), 
shows a peak at 6200 K for the HR compilation data while 
the MR measurements has a gradual increasing from 8000 
to 4200 K. There is a wide coverage over the whole MILES 
metallicity scale for both HR and MR data sets. In the pro- 
jection [Fe/H] vs. Toft, Fig. 7-(b), we notice a predomina- 
tion of the HR data between 5000 and 6000 K. In the plane 
log g vs. [Fe/H], Fig. 7-(c), the dwarfs are well covered by 
the HR data. The MR distribution dominates in the metal- 
rich regime and the HR data compilation dominates in the 
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Table 4. The MILES [Mg/Fe] catalogue for field stars. First column presents the star identification in the MILES database. The 
identification in the CaT library (Cenarro et al. 2002) is in the second column. The stars' names in other catalogues are shown in the 
third column. The stellar photospheric parameters in MILES are listed from the 4*'^ to 6"^ columns. [Mg/Fe] is shown in the seventh 
column together with its error in the eighth column. Notes about the source of each [Mg/Fc] measurement are written in the last column 

identifying its origin, i.e. from the high-resolution compilation (HR) or our mid-resolution measurements (mr), as well as the HR work, 
BM05 or other(s) listed in Table 1, the Mg feature(s) adopted in each stellar MR measurement, and if the mid-resolution measurement 
represents an a-enhancement model atmosphere extrapolation as described in Sect. 3.5 (designated by *). The full table is only available 
in electronic form. 



#MILES 


#CaT 


Star name 


Teff 


log 9 


[Fe/H] 


[Mg/Fe] 


CT[Mg/Fe] 


Notes 








(K) 




(dex) 


(dex) 


(dex) 




0081F 




BD-010306 


5650 


4.40 


-0.90 


-1-0.40 


0.05 


HR BM05 


0266F 




BD-011792 


4948 


3.05 


-1.05 


-1-0.71 


0.20 


HR T98 


0505F 


677 


BD+012916 


4238 


+0.35 


-1.47 


-1-0.37 


0.20 


HR T98 


0329F 




BD-032525 


5750 


3.60 


-1.90 


-1-0.41 


0.07 


HR BM05 


0777F 




BD+044551 


5770 


3.87 


-1.62 


+0.42 


0.07 


HR BM05 


0327F 




BD-052678 


5429 


4.43 


-2.14 


+0.41 


0.10 


mr BothMg 


0569F 




BD-f053080 


4832 


4.00 


-0.88 


+0.55 


0.10 


mr BothMg 


0142F 




BD-1-060648 


4400 


1.02 


-2.10 


+0.52 


0.20 


HR T98 


0144F 




BD-060855 


5283 


4.50 


-0.70 


-0.01 


0.05 


HR BM05 


0537F 




BD-l-062986 


4450 


4.80 


-0.30 


+0.06 


0.13 


mr Mg5183 



Table 5. The MILES [Mg/Fe] catalogue for star cluster stars (as in Table 4). The cluster names and types are in 3'''^ and 4*^ columns. 
The full table is only available in electronic form. 



#MILES 


#CaT 


Cluster name 


Type 


Star name 


Teff 


log 9 


[Fc/H] 


[Mg/Fe] 


(T[Mg/Fe] 


Notes 












(K) 




(dex) 


(dex) 


(dex) 




0920C 


006 


Coma Ber 


open 


HD107276 


7972 


4.21 


-0.05 


+0.29 


0.15 


mr Mg5528 


0921C 


007 


Coma Ber 


open 


HD107513 


7409 


4.25 


-0.05 


+0.32 


0.15 


mr Mg5528 


0901C 


016 


Hyades 


open 


HD025825 


5992 


4.41 


+0.13 


-0.09 


0.10 


mr BothMg 


0902C 


017 


Hyades 


open 


HD026736 


5657 


4.45 


+0.13 


-0.24 


0.10 


mr BothMg* 


0904C 




Hyades 


open 


HD027383 


6098 


4.28 


+0.13 


+0.03 


0.10 


mr BothMg 


0905C 


023 


Hyades 


open 


HD027524 


6622 


4.28 


+0.13 


-0.16 


0.15 


mr Mg5528* 


0906C 


025 


Hyades 


open 


HD027561 


6742 


4.24 


+0.13 


+0.03 


0.20 


HR T98 


0907C 




Hyades 


open 


HD027962 


8850 


3.80 


+0.13 


-0.09 


0.20 


HR T98 


0908C 


029 


Hyades 


open 


HD028483 


6486 


4.30 


+0.13 


+0.01 


0.20 


HR T98 


0909C 




Hyades 


open 


HD028546 


7626 


4.11 


+0.13 


-0.23 


0.15 


mr Mg5528* 



metal-poor regime, complementing each other well (see also 
Figs. 9-(c) and 11). Whilst the HR measurements provide 
mostly data for MILES dwarfs, reaching the completeness 
maximum around log g = 4.0 (Fig. 8-(b)), our MR mea^- 
surements contribute significantly to giants, with a gradual 
decrease from log g = 2.5 to 0.0. 

Specifically for the mid-resolution measurements that 
were done within the limits of MARCS 2008 grid as shown 
in Figs. 8 and 9, the highest completeness over the Teff scale 
occurs around 4000 K whilst the smallest is at 6000 K. The 
maximum coverage in the gravity scale occurs at log g = 2.5 
and the minimum is around log g = 4.0. The maximum of 
completeness over the metallicity scale is at —0.40 dex and 
the minimum occurs at [Fe/H] = —2.4 dex. 

Figure 9 presents the MILES catalogue's [Mg/Fe] dis- 
tribution itself, which is highly asymmetric around the solar 
ratio showing a sharp decline towards negative values and a 
shallow decreasing towards over-enhanced ratios (skewness 
= +0.129). The average of [Mg/Fe] is +0.16 dex having a 
standard deviation of 0.22 dex. The median of distribution is 
+0.12 dex. The distribution of our MR measurements match 



well the HR data distribution (both sets on a same homoge- 
neous scale). The median and average of HR data are +0.16 
dex and +0.19 dex (la = 0.19 dex) respectively, whilst they 
are +0.09 dex and +0.14 dex {la = 0.24 dex) for the MR 
measurements. The positive asymmetry also exists in both 
distributions: the HR data have skewness equals to +0.081 
and for the MR measurements it is +0.220. The difference 
between them is that the HR data present a peaked distri- 
bution (kurtosis = +0.178) whilst the MR values show a less 
peaked distribution (kurtosis = -0.279). The MILES stars, 
therefore, now have [Mg/Fe] measurements covering a range 
that is not restricted to the solar abundance ratio. 

Figure 10 plots [Mg/Fo] as function of [Fo/H]. Wo can 
affirm that the MR measurements follow very well the solar 
neighbourhood global pattern of the HR data, i.e. our de- 
terminations at medium spectral resolution statistically re- 
cover with acceptable accuracy this abundance ratio along 
the whole metallicity scale. At a given [Fc/H], the scatter 
of [Mg/Fe] is slightly larger when it is measured by using 
Mg5183 than when Mg5528 or both features combined are 
adopted. However, the systematic error on the [Mg/Fe] value 
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Figure 7. Coverage of the MILES stars with [Mg/Fe] in the 
photospheric parameter space: (a) modified H-R diagram log g 
vs. Teflf (top panel), (b) [Fe/H] vs. T^ff plane (middle panel), 
and (c) log g vs. [Fe/H] projection (bottom panel). The MILES 
stars with abundance ratios from HR studies are shown as black 
filled symbols, and the stars with MR, measurements from this 
work as colour filled symbols (blue designating determinations 
based on the Mg5183 feature only, red on Mg5528 only and 
green on both Mg features combined). The library stars with- 
out [Mg/Fe] are drawn as open symbols. The MARCS 2008 grid 
extension is also represented in each panel. 



Figure 8. Distribution of the MILES stars with [Mg/Fe] over 
the photospheric parameter space: (a) effective temperature scale 
(top panel), (b) surface gravity scale (middle panel), and (c) 
mctallicity cscale (bottom panel). The percentage number per 
bin of the MILES stars with abundance ratios from HR studies 
is shown by blue dashed line, and the percentage number of stars 
with our MR measurements by the red dotted line. The sum of 
both is represented by the black solid line. The adopted bins are 
equal to 4 times the parameter uncertainties. The coverage of the 
MARCS 2008 grid of model atmospheres is also drawn in each 
panel. 
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HR compilation 

our UR measurements 

sum 



(0 300 



^ 200 



2 100 



median = -HO. 12 
mean = +0.16 
Sigma = 0.22 
skewness = +0.129 
kurtosis = -0.110 



Table 6. Stellar parameter coverage of the MILES [Mg/Fe] cata- 
logue for the high-resolution and mid- resolution data distinguish- 
ing dwarfs (log g ^ 3.0) and giants (log g < 3.0). 



0.00 
[Mg/Fe] 



Class 


Toff 


log 9 


[Fe/H] 


[Mg/Fc] 




(K) 




(dox) 


(dox) 


HR 


Dwarfs 
Giants 


4342,11704 
3902,6666 


3.0,4.80 
0.0,2.98 


-2.53, +0.60 
-2.64,+0.10 


-0.54,+0.74 
-0.21, +0.65 


MR 


Dwarfs 
Giants 


3330,7972 
3600,7636 


3.0,5.08 
0.0,2.99 


-2.86,+0.41 
-2.50, +0.64 


-0.36, +0.73 
-0.47,+0.67 



Figure 9. Distribution of abundance ratios over the MILES 
[Mg/Fe] catalogue built up in this work. The histogram of stars 
with data from the HR compilation is shown by blue dashed line, 
and the distribution of our MR measurements by the red dotted 
line. The sum of both data sets is plotted by the black solid 
line. The adopted bin is two times the average uncertainty of 
[Mg/Fe] over the whole catalogue (2cr[Mg/Fe] = 0.21 dex). The 
median, mean, sigma, skewness and kurtosis of the accumulated 
distribution are also informed. 



Table 7. The distribution of HR, and MR measurements around 
dwarfs (log g > 3.0) and giants (log g < 3.0) in the MILES [Mg/Fe] 
catalogue. The weighted average uncertainty (T[Mg/Fe] is also pre- 
sented for each group of measurements. 



Data source Dwarfs Giants Total crfMg/Fc] 

dex 



HR 263 52 315 0.09 



is smaller when the Mg5183 feature is used (see also Sect. 
3.5). In the MILES data set there is a mix of stars from dif- 
ferent kinematic populations of our Galaxy distributed over 
the thin and thick discs as well as the halo. Several recent 
studies based on homogeneous HR spectroscopic analyses 
have shown the intrinsic dispersion of [Mg/Fe] or [a/Fe] at 
a fixed [Fe/H] inside a particular disc population seems to be 
really small indeed (e.g. Chen et al. 2000; Mishenina et al. 
2004; Bensby et al. 2005; Reddy, Lambert & AUende Prieto 
2006; Bensby et al. 2010; Nissen & Schuster 2010) but Neves 
et al. (2009) have found the opposite result. On the other 
hand, halo stars ([Fe/H] ^ —1.0 dex) exhibit great spread 
in the plane [Mg/Fe] vs. [Fe/H] (e.g. Stephens & Ann Mer- 
chant 2002; Borkova & Marsakov 2005). It is not the scope of 
the current work to explore the details about the elemental 
abundances over the Galaxy's kinematic components. These 
issues may be the central subject of a future work. 

Figure 11 shows that there is good coverage of [Mg/Fe] 
over the Tofr and log g scales (uniformly from 4000 to 5500 K 
and nearly uniform along whole gravity scale) , with poorest 
completeness at the lowest and highest temperature ranges. 
In addition there is a dearth of stars with sub-solar [Mg/Fe] 
around Teft = 6000 K, and for giant stars with low log g 
values. 

Table 6 shows the stellar parameter coverage of the 
MILES [Mg/Fe] catalogue over the HR and MR data, and 
dwarfs and giants as well. Table 7 summarizes the cata- 
logued data presenting the number of HR and MR measure- 
ments around dwarfs and giants together with their uncer- 
tainties. 



Mg5183 


23 


62 


85 


0.13 


Mg5528 


31 


69 


100 


0.15 


Both features 


96 


156 


252 


0.10 


MR sum 


150 


287 


437 


0.12 


Total sum 


411 


341 


752 


0.105 



5 COMPARISON WITH THEORETICAL 
MODEL PREDICTIONS 

In previous stellar population studies, there have been at- 
tempts to account for variations in element abundance ratios 
by the use of response functions that give the dependence 
of line-strength indices on just a single chemical element. 
These response functions are calculated by adopting model 
atmospheres of, usually, just 3 or 4 different stars. Examples 
include the models of Weiss, Peletier & Matteucci (1995), 
Tripicco & BeU (1995), Korn, Maraston & Thomas (2005) 
- hereafter K05, Coelho et al. (2005), and Serven, Worthey 

6 Briley (2005) {[a/Fe] = 0.0 and +0.3 models, effects on 
many spectral lines of many elements tested individually). 
Such models have been used in various studies to interpret 
the spectra of galaxies and globular clusters (e.g. Vazdekis 
et al. 1997, Trager et al. 2000b, Proctor & Sansom 2002, 
Denicolo et al. 2005, Lee & Worthey 2005, Schiavon 2007, 
Coelho et al. 2007, Pipino et al. 2009b, and Smith, Lucey 
& Hudson 2009). The accuracy of those response functions, 
however, have not been calibrated or tested empirically. 

In the present paper, variations in [a/Fe] are charac- 
terised by measuring Mg abundances in real stars and using 
this to represent a-elements in general. In this section, we 
compare the effects on Lick indices of these [a/Fe] ratios 
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Figure 10. [Mg/Fe] as a function of [Fe/H] sliowing all MILES stars with [Mg/Fe] from the HR compilation as black symbols and MR 
measurements as colour symbols (blue designating determinations based on the Mg5183 feature only, red on Mg5528 only and green on 
both Mg features combined). Dwarfs (log g 3.0) are shown as circles and giants as triangles. The weighted average uncertainties for 
each data group (HR, Mg5183, Mg5528 and both features combined) are illustrated on the bottom left corner. Three chemically peculiar 
stars are also identified. 



derived from observations with theoretical predictions for 
how Lick indices are expected to change with variations in 
[a/Fe]. The models of K05 are used for this comparison, 
since they were used in many of the above referenced stud- 
ies of galaxies. Equation 7 from Thomas, Maraston & Ben- 
der (2003) shows how we can predict changes in spectral 
line indices with changing composition for lines that tend 
to zero strength as the element abundance dominating that 
line tends to zero. Other line indices (e.g. H7, H^, G4300, 
Fe4383) take both negative and positive values, due to lack 
of a real continuum definition in complex star spectra. This 
particularly affects indices in the blue part of the spectrum 
where the continuum level changes rapidly with wavelength 
in long-lived, late-type stars. For these indices we adopt the 
formalism given in equation 3 of K05, which modifies fluxes 
rather than line strengths. For molecular bands and nega- 
tive going lines differences in indices are compared, whereas 
for positive absorption lines ratios are used for comparisons. 

The catalogue of MILES atmospheric parameters (Ce- 
narro et al. 2007) was cross-correlated with the measured 
[Mg/Fe] cases and with the models of K05 to find samples 
of stars useful for comparing observations with models, in 
order to test their agreement. There are then 31 stars in the 



MILES library whose surface temperature and gravity are 
the same as that of the turn-off star model (Toh = 6200 K, 
log (7 = 4.1) given in Table 13 of K05, within observational 
errors (ATgff = ±100 K, Alogg = ±0.2). Amongst these 
31 stars are two that also have the same chemical compo- 
sition as the model star ([Fe/H]=0, [a/Fe]=0), within ob- 
servational errors (A[Fe/II] = ±0.1 dex, A[Q/Fe] = ±0.06 
dex). This allows us to use one of these two stars as a base 
with which to normalise the other stars in order to show 
how changes in chemistry affect changes in indices, in a rel- 
ative way. Model stars can then be generated to match the 
31 MILES stars and normalised by the model given in Table 
13 of K05. To generate the model star indices the equations 
were first applied to correct to a specific overall metallicity 
(using column 14 of K05), then the equations were applied 
again to correct to a specific [a/Fe] ratio, modifying for all 
the a-elements modelled by K05. Similarly, 13 cool giant 
stars and 7 cool dwarfs can be compared using Tables 14 
and 12 from K05 respectively, for solar composition mod- 
els. In this way we can compare normalised observations 
with normalised model predictions to see if the observations 
agree with previously used methods of varying [a/Fe] ratios 
in stellar population studies. 
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Figure 11. Distribution of [Mg/Fc] in our catalogue over the 
scales of effective temperature (truncated at 9000 K) in the top 
panel (a), and surface gravity in the bottom panel (b). The 
notations are those adopted in Fig. 10. 



Examples of these comparisons are shown in Fig. 12 for 
Fe and Mg sensitive indices. Fe sensitive features in gen- 
eral behave as expected in that the observed changes agree 
well with the predicted ones. Mg sensitive features also show 
quite a good one-to-one agreement, but with scatter in ex- 
cess of that expected from the observational errors. There is 
some suggestion of a slight systematic deviation below the 
one-to-one line in the case if Mg b in turn-off and cool-dwarf 
stars. These deviations will be explored in future work. In 
general we see from Fig. 12 that the observed [Mg/Fe] abun- 
dances reported in this paper show the trends expected for 
Mg and Fe sensitive features, when compared to models from 
K05. Other indices and comparisons with models will be dis- 
cussed more extensively in a future paper. 



6 FURTHER APPLICATIONS TO THE 

ANALYSIS OF STELLAR POPULATIONS 

From now on, it will be possible to build up new semi- 
empirical simple stellar population (SSP) models by adopt- 



ing the MILES star spectrum library in order to more 
confidentially cover a range of values of metallicity and 
magnesium-to-iron abundance ratio for some ages greater 
than 1 up to 14 Gyrs. This is one of the important further ap- 
plications based on the results of this work that is scheduled 
by our group which has experience in studying stellar pop- 
ulations. The MILES stars can be now selectively collected 
from the library taking into account their characteristics in 
a more extensive parameter space, i.e. [Fe/H], [Mg/Fe], log 
g and Toff, to represent different SSPs for given sets of age, 
[Fe/H] and [a/Fe]. Magnesium may be considered a proxy 
of the a-elements and consequently [a/Fe] might be repre- 
sented by [Mg/Fe] . 

There are some caveats to be aware of for using the 
present catalogue results. The measurements made and com- 
piled in the present catalogue represent [Mg/Fe] well, as our 
tests of the spectral measurements show, considering all a- 
element or only Mg variations. However, it is important to 
be aware that not all a-elements may behave in exactly the 
same way in different populations. This needs to be con- 
sidered when applying the present catalogue to stellar pop- 
ulation studies. Other caveats are that [Mg/Fe] obtained 
through the MR calibrations applied for stars whose pa- 
rameters lay outside the control sample coverage might be 
more uncertain than the other determinations, and that we 
warn the abundance ratios derived from the a-enhancement 
model atmosphere extrapolations should be used with cer- 
tain precaution too. Caveats aside, making the approxima- 
tion that [Mg/Fe] can be used to represent [a/Fe] is a sig- 
nificant improvement over the scaled solar assumption only. 
There is a great deal of interest in uncovering the informa- 
tion contained in non-solar abundance ratios. Therefore we 
next illustrate how the catalogue may be used to generate 
SSPs with empirically determined non-solar [a/Fe] abun- 
dance ratios. 

In this further step of our work, it will be necessary 
to take into account reliable cross-matching between the- 
oretical isochrones for non-solar ratios and real stars. Ba- 
sically, the dwarf and giant stars with known Mg/Fe ratio 
that are collected from MILES in order to represent a given 
simple stellar population must be selected to precisely have 
on average the SSP's [Fe/H] and [Mg/Fe]. Moreover, the 
selected stars must be well sampled along the main evo- 
lutionary stages of an isochrone. The stars must be suffi- 
ciently close to the isochrones taking into account the errors 
in log g and Teff in order to be included and weighted in 
the computation of SSP integrated colours and spectra. It 
will be also necessary to generate additional stellar spectra 
for some non-completely-represented stages by interpolating 
the MILES spectra in its four-dimensional parameter space. 
Other approaches can be employed to extend the coverage 
of the SSP sets (age, [Fe/H], [a/Fe]) such as evaluating the 
behaviour of integrated observational properties of the SSPs 
at specified metallicities as a function of [a/Fe] to computing 
and applying corrections to the SSP observables for a larger 
range of parameters. Therefore we will be able to construct 
a large set of single-age single-metallicity single-a-enhanced 
stellar population models. This will open new prospects for 
SSP modelling and evolutionary population synthesis. 

Table 8 presents a set of age, [Fe/H] and [a/Fe] combi- 
nations (36 in number) for which there are sufficient number 
of dwarfs and giants in our MILES [Mg/Fe] catalogue to be 
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Figure 12. Normalised observed versus normalised model indices for different stars in the MILES library. The normalisation of the 
observations is achieved using MILES stars ■with the same photospheric parameters as for the models of K05 within errors. These 
normalising stars are labelled with 'base' in the plots. The top row shows MILES stars plotted against cool dwarf star models derived 
from Table 12 of K05; the central row shows MILES stars plotted against turn-off star models derived from Table 13 of K05 and the 
lower row shows MILES stars plotted against cool giant star models derived from Table 14 of K05. The models are normalised by the 
values given in tables 12 to 14 of K05. Models arc normalised in the same way as the data, as differences or ratios (see Sect. 5 for a 
description). Stars in a given plot all have the same effective temperature and surface gravity as that of the relevant model star, within 
the errors. The chemistry ([Fe/H] and [a/Fe]) is varied. 



more consistently selected to construct SSPs by adopting 
scaled solar and a-enhanced isochrones with those proper- 
ties. Basically, in this apijroacli, the selection of MILES stars 
for each metallicity value is done taking into account Icr or 
2a variation depending on the star's position in the H-R 
diagram. The stars are also collected having [Mg/Fe] val- 
ues around IcrfMg/Fe] and 2(7[Mg/Fe] of each isochrone's a- 
enhancement. The weighted uncertainty of the abundance 
ratio in the whole MILES [Mg/Fe] catalogue is 0.105 dex 
(Icr) and cr[Fe/H] = 0.10 dex in MILES. We also consider 
variations in the isochrones' ages, which were estimated ba- 
sically through plots of isochrones with distinct ages that 
are not shown in this section. 

The Dartmouth isochrone models (Dotter et al, 2008) 
have been adopted as reference to match the stars' posi- 
tions with isochrones in the HR diagram log g versus Tea. 
The Dartmouth models is a collection of scaled solar and a- 
enhanced isochrones that spans a range of [Fe/H] from —2.5 
to +0.5 dex, [q/Fc] from -0.2 to +0.8 dex (for [Fe/H] 
0.0 dex) or -0.2 to +0.2 dex (for [Fe/H] > 0.0 dex), with 
0.2 dex steps, and initial helium mass fractions from Y = 
0.245 to 0.40. Their stellar evolution tracks were computed 



for masses from 0.1 to 4 M©, allowing isochrones to be gen- 
erated for ages as young as 250 Myr up to as old as 15 Gyr. 

Figure 13 shows examples of isochrone-based plots by 
using Dartmouth (Dotter et al. 2008) and BaSTI (Pietrin- 
ferni et al. 2004) models for 4 Gyr, [Fe/H] around zero and 
three distinct a-enhancements. 

The BaSTI scaled solar database (Pietrinferni et al. 

2004) covers stellar evolution models for masses between 
0.5 and 10 M© in a wide metallicity range (10 values of 
[Fe/H] from -2.27 to +0.40 dex). The initial He mass frac- 
tion ranges from Y = 0.245, for the more metal-poor com- 
position, up to 0.303 for the more metal-rich one. For each 
adopted chemical composition, the evolutionary models were 
computed without (called canonical models) and with over- 
shooting from the Schwarzschild boundary of the convective 
cores during the central H-burning phase. The stellar mod- 
els are used to compute isochrones in a wide age range, from 
30 Myr up to 15 Gyr. The overshooting models provide a 
better match to the observations at [Fe/H] around solar, and 
for ages equal and higher than 4 Gyr. Besides these models, 
BaSTI a-enhanced models were computed for [a/Fe] fixed 
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at +0.4 dex and 11 values of iron metallicity between —2.62 
and +0.05 dex (Pietrinferni et al. 2006). 

Figure 14 presents a good example of cross-matching 
between a group of MILES stars and an isochrone to pro- 
ceed to a reliable semi-empirieal modelling of a single-age 
single-metallicity single-a-enhanced stellar population. The 
isochrone-based plot of this figure on the log g vs. T^s plane 
is done for 6 Gyr, [Fe/H] = -0.4 dex and [o/Fe] = +0.2 dex 
based on a Dartmouth theoretical model. The MILES stars 
have been carefully chosen with metallicity and [Mg/Fe] 
around the isochrone values within la. Furthermore, the re- 
strictive matching of the stars' positions with the isochrone 
has considered the MILES errors in log g and Tefi within 
la too. We have also computed the averages of [Fe/H] and 
[Mg/Fe] for the selected stars to check if there is agreement 
with the model's values. In this cross-matching SSP-MILES, 
9 dwarfs (log g ^ 3.7) and 21 giants have been selected 
whose averages [Fe/H] and [.Mg/Fo] arc -0.40 and +0.19 
dex respectively, in excellent accordance to the SSP's pa- 
rameters (the standard deviations are iT[Fe/H] = 0.07 dex 
and (T[Mg/Fe] = 0.06 dex). The stars are roughly well sam- 
pled along the isochrone although there are some empty 
places that might be filled through interpolations applied in 
the library parameter space. This procedure must make the 
stars' distribution more uniform around the model chemistry 
shown for instance in the plane [Fe/H] vs. [Mg/Fe], refining 
in this way the computation of the SSP observables. 

We have also begun a study of the dependence of ab- 
sorption line indices on [Mg/Fe], focusing on the observed 
behaviour of some indices of the Lick System as a function 
of the photospheric parameters, as described in Sect. 5. We 
intend to compute semi-empirical fitting functions for the 
main Lick indices in order to improve and extend the SSP 
models for different a-enhancements. 

The predictions of new semi-empirical SSP models will 
be compared with the observables of distinct composite stel- 
lar systems such as globular clusters, dwarf galaxies, ellip- 
ticals and spiral bulges. Consequently, the models will be 
very useful to understand their star formation histories and 
chemical evolutions. 

Moreover, we can test and apply the same approach 
of this work to obtaining the calcium abundances for the 
MILES stars. Indeed, there are interesting questions about 
how the calcium-enhancement behaves in several compos- 
ite stellar systems relative to other a-elements like magne- 
sium. For instance. Smith et al. (2009) found for a sample 
of 147 red-sequence galaxies from the Coma cluster and the 
Shapley Supercluster that the [Ca/Fe] ratio is positively cor- 
related with the velocity dispersion, at fixed [Fe/H], how- 
ever its dependence is significantly less steep than that of 
[Mg/Fe] . On the other hand, Pipino et al. (2009a) obtained 
that the [Ca/Fe]-mass relation is naturally explained by such 
a standard galactic chemical evolution model, and explained 
that the observed under-abundance of Ca with respect to Mg 
can be attributed to the different contributions from Type 
la and Type II supernovae to the nucleosynthesis of these 
two elements. 

Additional applications of the present [Mg/Fe] cata- 
logue will potentially improve areas of our understanding 
of stellar atmospheres and spectral flux distributions from 
different types of stars present in the MILES library. 



Table 8. Sets of age, metallicity [Fe/H] and [Mg/Fe] (repre- 
senting a-enhancement) for which there are reasonable number 
of dwarf and giant stars in MILES to build up semi-empirical 
SSP models. The star counts are shown for the main sequence 

(MS), main sequence turn-off (TO), sub-giants (SG) and red gi- 
ant branch (RGB) stages. The Dartmouth isochrones (Dotter et 
al. 2008) have been adopted as reference for this purpose. The 
stars have been selected in the [Fe/H]- [Mg/Fe] parameter space 
assuming specified ranges around these quantities, whose crite- 
rion of choice was such that the sum of stars in the MS + TO 
stages ^ 5 and in the SG + RGB stages ^ 10. The weighted 
average iT[Mg/Fe] is 0.105 dex in our catalogue and <7[Fe/H] is 
0.10 dex in MILES. Isochrone-based plots for the combinations 
assigned by asterisks are presented in Figs. 13 and 14. 
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Figure 13. Isochrone-based plots on the HR diagram log g vs. 
Tcff for age 4 Gyr, [Fe/H] = 0.0 dex and tiiree different [a/Fe] 
(-0.2, 0.0 and +0.2 dex), by adopting MILES stars witli [Fe/H] = 
0.0 dex and [Mg/Fe] around these o-enhancements with 1 and 2 
standard deviations (tr): top panel (a) for [Mg/Fc] = —0.2 dex, 
middle panel (b) for [Mg/Fe] = 0.0 dex, and bottom panel 
(c) for [Mg/Fe] = +0.2 dex. The Dartmouth isochrones (Dotter 
et al. 2008) are drawn as thick blue lines. In the panel (b), a 
BaSTI overshooting scaled solar model for [Fe/H] = +0.06 dex 
(Pietrinferni et al. 2006) is also drawn (cyan crosses). The library 
stars are shown as filled black squares for which [Mg/Fe] has Icr 
precision and as open red squares for 2(t. 
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Figure 14. Example of a cross-matching between a selected 
group of MILES stars and an isochrone for SSP modelling. Top 
panel: isochrone-based plot on the H-R diagram log g vs. T^ff 
for 6 Gyr, [Fe/H] = -0.4 dex and [a/Fe] = +0.2 dex by collect- 
ing MILES stars with metallicity and [Mg/Fe] around these val- 
ues with la deviation that restrictively match the correspondent 
Dartmouth isochrone (Dotter et al. 2008) taking into account the 
MILES errors in log g and Tj.jj within Icr too (error bars placed in 
the bottom left corner) . The dwarfs are represented as black filled 
circles and giants as red open triangles in all plots, where dwarfs 
have log g ^ 3.7, covering the stages MS, TO and SG listed in 
Table 8, and giants log g < 3.7. Bottom left panel: [Mg/Fc] 
vs. [Fe/H] plot for the selected dwarfs and giants, in which the 
two blue thick lines show the stars' average values of [Fe/H] and 
[Mg/Fe] (the line lengths arc equal to the standard deviations of 
the averages). Bottom right panel: [Mg/Fe] vs. T^ff plot for 
the selected dwarfs and giants, in which the horizontal blue thick 
line shows the stars' average [Mg/Fe] over the whole stellar sam- 
ple temperature scale. The averages of [Fe/H] and [Mg/Fe] are 
written on the top of left and right panels respectively. 



7 CONCLUSIONS AND SUMMARY 

We have obtained [Mg/Fe] for 76.3% of the MILES steUar 
library (411 dwarfs and 341 giants, 76% and 77% of the 
total respectively), suitable for SSP modelling. The weighted 
average uncertainty cr[Mg/Fe] is 0.105 dex over this MILES 
[Mg/Fe] catalogue. 

Compilation of high spectral resolution [Mg/Fe] abun- 
dance ratios in the literature was extremely useful in defin- 
ing a uniform scale for [Mg/Fe] and in obtaining an exten- 
sive reference sample for the calibration of abundance ratios 
measured in our work at medium resolution. We emphasize 
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that the cahbration of mid-resolution measurements is an 
important step to be done in the whole process to achieve 
reliable results in a homogeneous reference system. 

A robust spectroscopic analysis was carried out using 
the MILES mid-resolution spectra and LTE spectral synthe- 
sis of two Mg features. Two methods were applied through 
an automatic process: pseudo equivalent width and line pro- 
file fitting. 

The typical error of [Mg/Fe] from the collected and cal- 
ibrated high-resolution measurements is 0.09 dcx, and the 
uncertainties from our MR analysis range from 0.10 to 0.15 
dex, with a weighted average of 0.12 dex. Thus we show that 
the accuracy of our measurements from MR spectra is quite 
acceptable, but not better than those from HR analyses. It 
is possible to measure element abundances in many more 
stars with such accuracy at MR when a large control sam- 
ple from HR measurements is adopted. Hence this catalogue 
of [Mg/Fc] measurements will be useful for a range of appli- 
cations for stellar population modelling and understanding 
stellar spectra. 

The pattern of [Mg/Fe] versus [Fe/H] found for the stars 
in the MILES library is as expected for stars in the so- 
lar neighbourhood, ranging from sub-solar [Mg/Fc] for high 
metallicity stars to super-solar [Mg/Fe] for low metallicity 
stars, as shown in Fig. 10. 

[Mg/Fe] measurements approximately characterise the 
alpha-to-iron ratios in stars. Although not all types of stellar 
populations would be well sampled by the stars in this cat- 
alogue, applications exploring the effects of non-solar abun- 
dance ratios will be possible for such objects as globular clus- 
ters, spiral galaxies, various types of low luminosity galaxies 
and dwarf galaxies. The importance of this is: (i) abundance 
patterns in stellar populations hold clues to their histories 
and (ii) the accuracies of previously used characterisations 
of abundance patterns, based on theoretical models, can now 
be tested. These applications will be followed up in future 
work. 

We also plan to use the MILES [Mg/Fe] catalogue to 
compute empirical and self consistent models of stellar pop- 
ulations with Mg/Fe different from solar for certain values of 
age and metallicity. These models will serve as a benchmark 
for other models based on theoretical libraries as they will 
allow calibration of the effects of uncertainties in the final 
predictions due to uncertainties in specific groups of stars. 
We will study empirically the dependences of Lick System 
indices as a function of the Mg/Fc ratio by adopting MILES 
stars with similar photosphcric parameters but showing dis- 
tinct [Mg/Fe] in order to help computing robust stellar fit- 
ting functions of line-strengths. Comparisons of empirical 
and theoretical line strengths will also be made (Sansom et 
al., in preparation). 
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Figure Al. Comparison of our HR data with the [Mg/Fe] com- 
pilation for the CaT library (Cenarro et al. 2009) that both have 
been independently calibrated to that scale of BM05 (whose mea- 
surements are represented by the filled black squares). The filled 
colour squares represent different data sets; blue for CGSOO, green 
for T98, yellow for FOO, and red for the stars with duplicated 
sources. The linear Isq fit [Fe/HJconOQ = A + B [Fc/HJmiles 
(solid blue line) is shown with parallel dashed blue lines illus- 
trating the 3(7 data clipping (its A and B, rms and correlation 
coefficient r are presented). The statistically representative in- 
verse expression [Fe/H]MILES = —A/B + 1/B [Fe/H]co09, based 
on a 95% t-test, is also shown on the top. The excluded data 
is represented by open colour squares (3-0" criterion and outliers 
from T98). 



APPENDIX B: 

In this appendix we investigate possible systematic differ- 
ences between the atmospheric parameters of the MILES 
catalogue (Cenarro et al. 2007), that have been adopted all 
over in the current work, and those of the BM05 sample. 
Just concerning the metallicity scale, comparisons are also 
made with the data of those consuhed HR works (see Sect. 
2). 

Figure Bl-(a)-(l) presents [Fe/H] from 12 consulted 
works compared with the MILES [Fe/H] scale. Lines of least- 
square {Isq) linear fittings [Fe/H]work = A + B [Fe/H]MiLES 
are shown in these plots. The adopted Isq method takes into 
account the errors in both variables by minimizing the sum 
of distances of all points to the line. Statistically representa- 
tive inverse linear transformation expressions [Fe/H] miles ~ 
—A/B + 1/B [Fe/H]work are obtained after applying a 95% 
confidence level t-test for each fit parameter {A^07 and/or 
B^ll). The fitted straight lines found are very close to the 
1:1 relationship. There are tiny systematic differences be- 
tween the MILES [Fe/H] scale and the scales from Fulbright 
(2000), Bensby et al. (2005), and Luck & Heiter (2005). How- 
ever, these differences are dominated by a few outliers, as 
can be seen in the panels (c) , (e) and (h) , that are specifically 
localized either at the metal-poor regime on the comparison 
with the Fulbright (2000) data, or over metal-rich stars on 
the comparisons with the samples of Bensby et al. (2005), 
and Luck & Heiter (2005). If we applied metallicity cor- 
rections for these samples, they would be smaller than the 
involved uncertainties, even for the works of Bensby et al. 
(2005), and Luck & Heiter (2005). Therefore no correction 
has been applied to the [Fe/H] from the consulted works be- 
cause no significant systematic deviations from the MILES 
[Fe/H] scale has been detected. 



APPENDIX A: 

This appendix shows the comparison of our [Mg/Fe] refer- 
ence scale with the CaT catalogue (Cenarro et al. 2009). 

Figure Al presents a direct comparison of the [Mg/Fe] 
high-resolution measurement compilation done by Cenarro 
et al. (2009) for their CaT library stars with the HR part of 
MILES [Mg/Fe] catalogue. A linear Isq fitting between their 
results and ours (also considering the errors in both variables 
and minimizing the distance along both directions) gives 
[Mg/Fe]caT = [Mg/Fe]ce09 = 0.00 -f 0.95 [Mg/Fe]MiLES 
with a data spreading of 0.03 dex {rms). This fit shows 
deviations, which are 0.035 dex at maximum, comfortably 
within the correspondent uncertainties of [Mg/Fe] in CaT 
(typically 0.14 dex) and MILES (0.09 dex on average). Two 
pairs of data were excluded by the 3-sigma clipping treat- 
ment taking into account their errors as well. Few outliers 
with high and low [Mg/Fe] values (all from T98) were also 
excluded from this comparison. The t-test with a 95% con- 
fidence level appoints that the relationship [Mg/Fe] miles 
versus [Mg/Fe] Ceog is not far from the 1:1 relation (see Fig. 
Al). Therefore there is a quite good agreement between both 
scales. 



For the comparison of BM05 and MILES TcS scales, the 
rms of 83 K for a linear relationship between the scales and 
the respective systematic deviation are comparable to the 
typical temperature uncertainty in MILES (lu = 100 K). In 
the range 4500-7000 K, for instance, the maximum absolute 
difference reaches around 55 K only. For the BM05-MILES 
log g scale comparison, the absolute difference in the interval 
3.0-5.0 gets a maximum of 0.13, which is smaller than the 
MILES log g uncertainty (la — 0.20). The rms of the linear 
relationship between the gravity scales is smaller than the 
gravity's uncertainty too (0.145). The angular coefficients 
of the linear Isq fits Tcff(BM05) vs. TeH (MILES) and log 
g(BM05) vs. log 5i(MILES) are different, respectively, 1% 
and 11% only from that of the 1:1 relation. Therefore the 
scales of TcH and log g in BM05 and MILES agree quite well 
between each other, and no correction has been applied to 
these stellar parameters too. 
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Figure Bl. [Fe/H]„ojij vs. [Fs/HJmiles- 12 panels (from a to 1) showing comparisons between the metaUicity scales of the consulted 
HR works (designation at the bottom of each panel) with the MILES one. The linear Isq fittings [Fe/H]„ork = A + B [Fe/H]MiLES with 
a 3-sigma data clipping are presented and illustrated by respectively solid blue line plus parallel dashed blue lines with excluded points 
in red. The constants A and B, rms and correlation coefficient r are shown on top of each plot as well as the statistically representative 
inverse expressions [Fc/HJmiles = —A/B + l/B [Fe/H]„orki based on 95% t-tests. 
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APPENDIX C: 

In this appendix, comparisons with stars cluster data from 
high-resolution studies are presented. 

There are 89 cluster stars in the MILES database (from 
9 open clusters and 8 globular clusters). Our mid-resolution 
measurements cover 65 cluster members (73% of them) in- 
cluding 16 out of 17 clusters of MILES; see Table 5 (Sect. 
4). Specifically for three star clusters presented in MILES 
(Hyadcs, M71 and NGC7789), wc have obtained [Mg/Fc] for 
a reasonable number of members from our spectral synthesis 
at mid-resolution (^10 stars). Average values arc presented 
in Table CI together with the results from other clusters for 
which we have done MR measurements for a minimum of 
two stars, providing an interesting quality test for our work. 
Table also shows average cluster values of [Mg/Fe] collected 
from HR studies, but their abundance ratios are not in the 
same uniform scale adopted in the current work. The litersr 
ture cluster averages [Mg/Fe] are computed from published 
ratios to be representative to the star sample analysed in 
our work for each cluster in terms of spectral type and lu- 
minosity class. 

For the open cluster Hyadcs only ([Fc/H] = +0.13 dex), 
we could compare our [Mg/Fe] MR measurements directly 
with the compiled HR data in a same uniform scale. The av- 
erage [Mg/Fe] computed from the MR measurements based 
on the two Mg features and including 10 members (dwarfs 
only with 5256 < Tea 7634 K) is -0.035 dex with a stan- 
dard deviation a = 0.169 dex, whilst the HR data have a 
average value —0.016 (Icr = 0.063 dex) from 3 dwarfs only 
(all from T98 with Teff = 6486, 6742 and 8850 K). The av- 
erage of [Mg/Fc] from both MR and HR data for Hyades is 
-0.030 (cr = 0.149 dex). Schuler, King & Lih-Sin (2009) has 
recently measured, by analysing high-S/N high-resolution 
spectra, [Mg/H] = -1-0.10 (cr = 0.02 dex) in 3 main sequence 
stars (Toff around 5600 K), getting [Mg/Fe] = -0.03 dex. A 
sample of 55 F-K dwarfs (with 4900 5? Tcff 6450 K) were 
spectroscopically analysed at HR and high-S/N by Paulson, 
Sneden & Cochran (2003). They obtained [Fe/H] = -1-0.13 ± 
0.01 dex and differentially [Mg/Fe] = -0.03 ± 0.04 dex. The 
spectroscopic analysis (based on HR and high-S/N spectra) 
of Yong, Lambert & AUcndc Pricto (2004) measures for 34 
Hyadcs dwarfs with 4700 < Teff ^ 6200 K cluster averages 
of [Fe/H] = -0.16 {a = 0.10 dex) and [Mg/Fe] « -0.09 dex. 
All these results corroborate our work. 

For the globular cluster M71 ([Fe/H] = -0.84 dex), in 
which there is a known internal spread of a/Fe ratio over 
the stars (Ramirez & Cohen 2002), we obtained an aver- 
age [Mg/Fe] abundance ratio equals to +0.485 (a = 0.104 
dex) based on measurements from both Mg features in 19 
cluster members, which arc all giants with 4014 ^ Toff ^ 
5123 K. [Mg/Fe] varies in M71 as Ramirez & Cohen (2002) 
measured in 24 giants that were spectroscopically analysed 
at HR. This star sample exhibits an average [Mg/Fe] equals 
to +0.36 (a — 0.09 dex) and a maximum internal spread 
of 0.18 dex. Melendez & Cohen (2009) recently proved, as 
concluded in other previous studies, the existence of two 
stellar populations in M71; one CN-weak with normal O, 
Na, Mg, and Al abundances plus a low isotope abundance 
ratio of ^"Mg/'^^Mg, and other CN-strong with enhanced Na 
and Al accompanied by lower O together with a higher ra^ 
tio ^®Mg/^"'Mg. However, they measured a small spread for 



[Mg/Fe] over 9 giants (0.10 dex at most) belonged to both in- 
ternal populations exhibiting an average around +0.20 dex. 

For the galactic cluster NGC7789 ([Fe/H] = -0.13 dex), 
the average [Mg/Fe] is +0.077 (cr = 0.092 dex) from our 
Mg5183/Mg5528 measurements in 13 giant stars (with 4020 
< Toff < 4952 K). A very recent work (Pancino et al. 2010) 
has measured for it an average [Fe/H] = +0.04 ± 0.07 dex 
with 1(7 dispersion of 0.10 dex and [Mg/Fe] = +0.22 ± 0.07 
dex (cr = 0.10 dex) based on spectroscopic measurements of 
3 red clump stars with high-quality spectra at HR, whilst 
previous studies at lower spectral resolution and through 
photometry-based techniques obtained [Fe/H] ~ —0.2 with 
[X/Fe] around zero for many elements (e.g. Priel et al. 2002, 
Pilachowski 1985) 

Concerning our star cluster data, the major conclusions 
are: (i) the standard deviations of computed average clus- 
ter values of [Mg/Fe] are comparable with the systematic 
uncertainties of our individual MR measurements as well as 
with the dispersion of the cluster average abundance ratios 
collected from the HR studies, and (ii) the [Mg/Fe] averages 
from the current work are in good agreement with measure- 
ments carried out with high-S/N high-resolution spectra in 
recent studies. 

This paper has been typeset from a T^^X/ IM^jX file prepared 
by the author. 
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Table CI. Averages values of [Mg/Fe] for star clusters in the MILES library computed from two or more star individual measurements 
at mid-resolution (mr) done in this work, [Mg/Fe]mr and standard deviation <T[Mg/Fe]mr, and respective average from the high-resolution 
(HR) studies, [Mg/Fe]HR and standard deviation (T[Mg/Fe]HR- The number of stars for each average (mr and HR) is informed in the 
columns six and ten respectively. The clusters' metallicities adopted in MILES are listed in the third column. The Mg features adopted 

in our MR measurements arc presented in the seventh column and the references consulted for the HR data arc lied in the last column 
(GMR08 for Gcbran, Monicr & Richard 2008, SKL09 for Schuler, King & Lih-Sin 2009, YLA04 for Yong, Lambert & AUcndc Pricto 
2004, PSC03 for Paulson, Sncdcn & Cochran 2003, T98 for Thcvenin 1998, KeOO for King et al. 2000, CM05 for Cohen & Melcndez 2005, 
Je05 for Johnson et al. 2005, Sc04 for Sneden et al. 2004, 896 for Shetrone 1996, MC09 for Melendez & Cohen 2009, RC02 for Ramirez 
& Cohen 2002, G089 for Gratton & Ortolani 1989, and PclO for Pancino et al. 2010, on the order they appear in the table). The HR 
data from T98 (Hyades) comes directly from the MILES [Mg/Fe] catalogue. Note: the others [Mg/Fe]HR are not transformed onto the 
catalogue's uniform scale. 



Cluster 


Type 


[Fe/H] 








Mg feature (s) 






Nhr 


Ref 


[Mg/Fe] 


<T [Mg/Fe] 


[Mg/Fe] HR 


(T [Mg/Fe] HR 






(dex) 


(dex) 


(dox) 






(dex) 


(dex) 






Coma Ber 


open 


-0.05 


-1-0.304 


0.022 


2 


Mg5528 


-1-0.26 


0.008 


2 


GMR08 


Hyades 


open 


-f0.13 


-0.035 


0.169 


10 


Both 


-0.03 


0.02 


3 


SKL09 
















-0.09 


0.10 


34 


YLA04 
















-0.03 


0.04 


55 


PSC03 
















-0.016 


0.063 


3 


T98 


Pleiades 


open 


-0.03 


-0.122 


0.083 


2 


Both 


-0.01 


0.06 


2 


KeOO 


M3 


globular 


-1.34 


-fO.299 


0.183 


3 


Both 


-1-0.41 


0.12 


13 


CM05 
















-1-0.17 


0.15 


77 


Je05 
















-1-0.22 


0.15 


23 


Se04 


M5 


globular 


-1.11 


-1-0.426 


0.169 


2 


Both 


-1-0.16 


0.07 


6 


S96 


M71 


globular 


-0.84 


-1-0.485 


0.104 


19 


Both 


-1-0.20 


0.10 


9 


MC09 
















-1-0.36 


0.09 


24 


RC02 
















-1-0.34 


0.08 


8 


S96 


M79 


globular 


-1.37 


-fO.493 


0.009 


2 


Mg5528 


-t-0.47 


0.45 


2 


G089 


M92 


globular 


-2.16 


-fO.485 


0.173 


2 


Mg5528 


-1-0.19 


0.19 


6 


S96 


NGC7789 


open 


-0.13 


-1-0.077 


0.092 


13 


Both 


-fO.22 


0.10 


3 


PelO 



